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Abstract 

In this paper, we investigate the effect of state and 
measurement variations on linear time-invariant state 
(LTI) observers’ performance robustness and error 
dynamics. These observers are widely used in control 
engineering, and, more recently, modern automotive 
powertrain applications. Their approach is analogous to 
Kalman Filters except that the perturbations are 
deterministic as opposed to stochastic. They have a 
proven and well documented degree of reliability and 
efficiency in the approximation of systems’ states. When 
well designed, they exhibit negligible steady state errors. 
In this study, we apply observer theory to an internal 
combustion engine (ICE) model. We also investigated 
state and sensor measurement disturbances that occur 
during common operation of the system. An analysis of 
the change in the observer error dynamics, and 
subsequent performance robustness is presented. 

Keywords: Linear observers; Control theory; Internal 
combustion engines modeling; State disturbance; Sensor 
disturbance; Estimation methods. 

1. Introduction 

There has been much work done on LTI state observers 
and there results are well known in control engineering. 
The concept was first introduced by Luenberger in 1971 
[1]. Observers were presented both for the continuous 
and discrete-time cases. Since then, state observers have 
revolutionized control theory by enabling the 
approximation of the state vector using readily available 
measurements of the system’s outputs. These outputs 
come directly from measurements. In some systems 
however, not all of the state variables can be easily 
observed through direct measurement. In such cases, 
costly measurement devices might also be needed to 
measure these states. In other systems, some of the state 
variables may be inaccessible for direct measurement. 
For the aforementioned cases, a state observer provides a 
good solution by approximating the states for which 
information would otherwise not be available. This 
makes the state observer inherently invaluable as a 



model-based method for state-feedback control of any 
control system. In the four decades following the 
introduction of the LTI state observer methodology, its 
use has been extended through various disciplines of 
engineering, including computer-aided software design 
and has been utilized in modem control stmctures to 
implement state observers in the manner depicted in Fig. 
1 . 

One of the recent fields that has been a beneficiary of the 
state observer is automotive powertrain. Specifically, we 
discuss its impact on ICE powertrain control. It is widely 
known that due to government mandates modem ICEs 
will have to improve their fuel efficiency and decrease 
their emissions if they are to remain the primary 
propulsion source [2]. As a result, new technologies, 
such as ‘engine start-stop’, and more modern control 
strategies, such as advanced idle-speed control schemes 
(ISC), have been employed [3-4]. The more popular 
method to solving efficiency issues has been to utilize the 
wealth of knowledge and proven effectiveness of 
model-based methods, which use state 
observers [5-11]. 




Figure 1: State observer feedback loop flow diagram 



The long established method of model-based controllers 
in the automotive industry to achieve fuel savings 
through better idle speed control (ISC) was first 
pioneered by Hrovate [12] in the mid 1990s. Around the 
same time, Zohdy and Karam began using optimal 
control to develop a more robust controller that addresses 
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the influences of the unmolded dynamics of ICS [13]. 
These unmolded dynamics include measurement noise, 
external disturbances, transient loads and model 
inaccuracies. To address these, the control community 
has developed methods such as Adaptive Controllers 
[14], nonlinear adaptive Fuzzy Logic [15], and Kalman 
filter [16] controllers. Kalman filters are most often used 
due to their inherent ability to deal with noise [17-18]. 
They have shown great success in many applications [19- 
23]. However, they are limited in many areas where the 
assumption that the noise has Gaussian distribution does 
not hold. Also, Kalman filters can be computationally 
complex to implement. Adaptive control schemes are 
limited by their need for high-fidelity models running in 
multi kHz switching frequency range and by the fact that 
they are highly nonlinear. These limitations have kept 
them out of the automotive market. 

One solution is the lower order, robust, linear time- 
invariant state observer model addressed in this paper 
which is computationally efficient enough to run in real- 
time in typical powertrain control modules (PCM). State 
observers are attractive in control schemes because they 
enable linear control schemes to be implemented on 
nonlinear systems. Unlike nonlinear control schemes, 
linear control schemes do not demand powerful 
processors; therefore, they minimize computational 
delays. While observer theory is well known in the 
literature, observers performance under non-ideal 
conditions and disturbance hasn’t been fully investigated. 
In this work, we include the influence of variations due to 
inaccuracies in the state space model and include 
measurement errors due to sensor errors. Model state 
disturbances and sensor measurement disturbances are 
investigated by analyzing their effects on observer error 
dynamics and overall state observer performance. 
Observer model designed to estimate the system states 
will be discussed in Section 2 and specifically applied to 
an internal combustion engine (ICE) model as introduced 
in [24] in Section 3. The results and explicit analysis of 
the findings are provided in Section 4 with our 
conclusion in Section 5. 

2. Observer Model 

In this section, a thorough derivation of the LTI state 
observer error dynamics is presented. The state observer 
is designed in order to reconstruct the states of a linear 
system [25-26]. The state space representation of the 
system includes the state and output equations, 
respectively presented in (1-2), 

x(t) = Ax(t) + Bu(t) + d (1) 

yif) - Cx(t) + rj (2) 

where X is the state vector of dimensions nxn and y is 
the output vector of dimension mxfl. The number of 
states in the state vector and output are denoted by n and 
m, respectively. The constant state, input and output 
matrices of appropriate dimensions are denoted by A, B 
and C , respectively. The inputs has dimensions rx\ , 
where V is the number of inputs to the system. This state 
space system includes the effect of the state disturbance, 



d , and sensor (output) noise, 77 . This is very similar to 

the state space model a Kalman filter would be applied to 
[27-28]. The effect of variations is of interest in this 
work and is represented in the model. The state space 
model mathematically demonstrates that dynamic 
changes in the state vector are affected by its current 
values at some time, t , the system’s input and 
disturbances on the state vector itself. The outputs of the 
system are measurements from the state vector, which 
take the errors of the measuring device into account. 

The state observer model is very similar to the 
state space model presented in (1-2) and is shown in (3-4). 

x(t) = Ax(t) + Bu(t) + L(y(t)-y(t )) (3) 

y(t) = Cx(t) (4) 

where x(t) denotes the estimate of the state vector and 
y(t) is the estimated output of the state observer model, 
which is comprised of the estimated state vector. It is 
easy to see that the observer state equation in (3) differs 
from (1) because the observer gain, L , is used as a 
correction term for the output error (y(t) — y(t)). One 
can notice that if the output error y(t) — y(t) is very 
small, then (3) and (1) are identical. However, if the 
output error is large, then L is used to correct the error in 
a way that will drive the overall state observer error 
dynamics to zero. 

There are well known methods to determine the observer 
gain, most notably the pole placement method [29-30] 
and the Linear Quadratic Regulator (LQR) method [31]. 
A block diagram of the state observer system is shown in 
Fig. 2. The state observer performance is determined by 
the error dynamics of the observer. A well designed 
observer has error dynamics that converge to zero 
quickly. How quickly depends on the criteria of the 
observer application. 



Figure 2: Block Diagram of state observer system 




The state observer error dynamics are derived as follows, 



e(t) = x(t) - x(t) 


(5) 


e(t) = x(t)-x(t) 


(6) 


e (t) = (A x (t) + B u (t) + d ) - (A x (t) + B u (t) 
+L (C x (t) + rj -C x (t))) 


(7) 


e (t) = A x (t) + B u (t) + d — Ax (t)-Bu(t) 
—LC x (t) — Lrj +LC x (t) 


(8) 


e(t) = (A- LC)(x(t) - x(t)) + d-Lri 


(9) 


e(t) -(A — LC) e(t) + d - Lri 


(10) 




ACSE Journal, ISSN: Print 1687-481 1 , Online 1687-482X, CD-ROM 1687-4838, ICGST LLC, Delaware, USA, December, 2014 



where (10) is the LTI state observer error dynamics. Our 
aim is to look at the effect variations of d and 77 have on 
the observer error dynamics and overall observer 
performance. 

3. Internal Combustion Engine 

The model used for the state observer error dynamics 
analysis is similar to the one presented by C.E. 
Baumgartner et al. in [24], where it was created and first 
implemented for ISC of the Dodge Neon (ICE). The 
model represents the intake dynamics of the engine in 
order to use engine speed to control idling mode. The 
initial model had three inputs and three states. Later, the 
model was modified by Omekanda and Zohdy in their 
work published to this publication [32]. It was reduced 
to only one input for simplicity and used for the design of 
a reduced order observer. The latter model is used here, 
with the input u(t ) and the state vector 
x(t) = [x x (t) x 2 (t ) x (t)] seen in Table 1 . 



Table 1 : Input and State vector 



u(t) : Commanded Throttle Angle (degrees) 



*1 (0 



Throttle angle (degrees) 



often too computationally complex to implement within 
the limitations of an automobile PCM. Our approach to 
the EIV problem is to design a computationally efficient 
controller robust enough to operate with error variations. 
Thus, we reduce the need to make assumptions about 
noise and make the controller more implementable in a 
modem PCM. In order to look into the effect of 
disturbance, d is assumed to be a random disturbance on 
the state that could be caused by faults or other 
unexpected variations. The output of the model is engine 
speed, and we depict sensor measurement disturbance as 
7 7- 




x 2 (t) : Intake manifold Pressure (Torr) 

x^it) : Engine speed (RPM) 



Figure 3: State Observer approximation of the ICE system with no 
disturbance 
Full Observer Error Dynamics 



It has been modified in such a way that an open loop 
throttle request of 7 degrees commands an idle speed of 
680 RPM. The state space of the ICE is presented as, 





-35 


0 


0 1 




35 


x(t) = 


4.212 


-3.822 


-6.552 


X(t) + 


0 




0 


3.2 


1.84 




0 



y(t) = [0 0 l]x(t) + »7 (12) 




An LTI state observer was designed for this ICE model 
when no state and measurement noise are taken into 
account. The gain, L, for this system was determined 
using the pole placement method. Figure 3 depicts the 
observer state approximation of the system states. Figure 
4 shows the state observer dynamics. The state observer 
has great dynamics, which converge very quickly and 
have no steady state error. During most of the first 
second of simulation, it appears that the intake manifold 
pressure is positive. Knowing this is physically 
impossible for a naturally aspirated ICE, we attribute this 
to a reduced order non bounded mathematical modeling 
error, resulting from an extremely transient step input of 
the original system of equations. 

To determine the model’s robustness and to evaluate the 
performance and effect on error dynamics, this paper will 
evaluate the ICE model’s response to measurement and 
modeling disturbance, often called error-in-variables 
(EIV). Much has been written on the difficulty of 
designing controllers to address EIV [33-35]. This leads 
to the belief that prior knowledge of states and several 
assumptions about the system must be made to control a 
system laden with EIVs [36-37]. These solutions are 



Figure 4: State observer error dynamics 



4. Results and Analysis 

The state disturbance is quantified as an assumed random 
disturbance bounded by a percent error (+/-) of the actual 
state (throttle angle, manifold pressure and engine speed), 
while the measurement disturbance is bounded by a 
percent error (+/-) of the engine speed sensor 
measurement. We first examined the effect of the state 
disturbance, d , and the measurement disturbance, 77 , by 
themselves, on the observer error dynamics and state 
observer performance. Then, we mixed different 
amounts of state and measurement disturbances and 
analyzed the resulting effects. The disturbances used for 
the singular disturbance tests and their results are shown 
in Table 2 below. 
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Table 2: Tests points for Figure 5 (a-d) 



Measurement Disturbance ( 77 )% 


State Disturbance (d) % 


0 


2 


0 


2 


0 


10 


0 


15 


2 


0 


5 


0 


10 


0 


15 


0 



Figures 5 (a-d) shows state observer error dynamics 
behavior caused by varying the state disturbance on the 
system by +/- 2, 5, 10 and 15%, respectively. There is no 
measurement disturbance assumed present in these tests. 
This means that engine speed measurements from the 
plant (ICE) used by the state observer to estimate the 
state vector are not corrupted. One should notice right 
away that the intake manifold pressure and engine speed 
have been scaled by a factor of 100. This was done to 
more easily demonstrate the instantaneous state and 
measurement disturbance included in the plots. It is 
immediately visible that the state observer is very 
sensitive to state disturbance. This is expected because 
an observer is a model-based method, which relies on an 
accurate model of the system to approximate the state 
vector. To account for state disturbance, we could adjust 
the gain in order to correct the error dynamics. Such 
methods have been implemented in literature by using 
gain scheduling schemes [38-45]. But in this study, we 
chose to leave the gains constant to illustrate the effects 
of the noise on the system response. As the state 
disturbance increased, the state observer error was unable 
to converge, and the oscillations and instability of the 
dynamics increased. In figure. 5a, the states try to 
converge but the observer seems to be marginally 
working. This means that it keeps the error within a 
bounded error range but never converges to zero or lost 
control. The throttle angle is the least affected by the 2% 
state disturbance. However, by figure. 5d the observer 
loses all controllability and the oscillations are as high as 
-10 degrees, -2000 Torr and -2000 RPM for the throttle 
angle, manifold pressure and engine speed, respectively. 



Error Dynamic Behavior(Throttle Angle) 
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| State Observer (Deg.) 


r~ " — 
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Figure 5b: 5% state disturbance 
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Figure 5a: 2% state disturbance 



In the next set of tests, no state disturbance was input into 
the system. Rather, the effect of measurement 
disturbance was analyzed. Our measurement disturbance 
was applied to the speed sensor, which is our output 
signal. Figure 6 (a-d) depict the error dynamics 
performances of each test when the measurement error on 
the speed was varied from +/- 2, 5, 10 and 15%, 
respectively. Compared to the effect of state disturbance 
only, the measurement disturbance seems to be less 
sensitive. This is evident in the lesser bounded 
oscillations at each disturbance percentage in comparison 
to those of Figure 5 (a-d). Here, the throttle angle 
approximation again seems to be the least sensitive to the 
disturbance. 

In Figure. 6 a, the error converges and only oscillates 
about zero during certain ranges of operation. Even in 
Figure. 6 b, small oscillations are evident and start to have 
larger bounds on the oscillations. From Figure. 6 c to 
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Figure. 6d, the error bounds for the oscillations barely 
demonstrate change. With greater than 5% measurement 
disturbance, the error dynamics become uncontrollable 
bounded by a certain error bound. 

Error Dynamic Behavior(Throttle Angle) 




Time (sec) 

Error Dynamic Behavior(Engine Speed) 




Figure 6a: 2% measurement disturbance 
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Figure 6b: 5% measurement disturbance 
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Figure 6c: 10% measurement disturbance 
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Figure 6d: 15% measurement disturbance 

A final set of conditions was when both measurement 
disturbance and state disturbance were introduced to the 



system together. Table three shows 4 points that where 
investigated. 



Table 3: Tests points for Figure 6 (a-d) 



Measurement Disturbance( 7 } ) % 


State Disturbance (d) % 


2 


2 


5 


5 


10 


10 


15 


15 



Through casual observations, one can see that these 
results are very similar to the state only classification of 
disturbances. These results confirm earlier findings that 
the system is much less sensitive, to an order of 10 
magnitudes, to measurement noise than state disturbance. 
Further analysis of these results was achieved by 
repeating the test with different sets of random noise. By 
capturing and comparing the state error between tests, it 
was shown that the relationship of the two statistically 
independent disturbances has a causal effect of the state 
observer. The variance of the state observer error was 
greater that with state disturbances alone, indicating that 
the two disturbances interacted to randomly increase or 
decrease their combined effect on the state error. The 
results of one of these random simulations are shown in 
figure 7. 
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Figure 7a: 2% state and measurement disturbance 

Error Dynamic BehaviorfThrottle Angle) 



| Stale Obseiver (Deg~)~| 






Time (sec) 

Error Dynamic Behavior( Intake Manifold Pressure) 



- 


1 


1 1 1 


State Obseiver (ToirxIOO) | 


- 


I 


I i i 


1 


) 1 2 3 4 5 

Time (sec) 

Error Dynamic Behavior(Engine Speed) 


- 


i 


i i i 


| State Obseiver (RPM xlOO) | 



Time (sec) 



Figure 7b: 5% state and measurement disturbance 
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Error Dynamic Behavior(Throttle Angle) 







Figure 7c: 10% state and measurement disturbance 




Time (sec) 



Figure 7d: 15% state and measurement disturbance 

5. Conclusion 

Our results demonstrate that state observers are very 
sensitive and susceptible to disturbance. This is in large 
part due to the fact that, unlike the Kalman filter, state 
observers do not take into account possible disturbances 
in their algorithm. On the contrary, any disturbance in 
the system directly affects the observer. However, state 
observers, as we demonstrated, are great state 
approximation tools for any system. In practice, they 
continue to be utilized because the control of systems is 
not based solely on the observed state(s) but also on the 
other sensor measurements and state observer 
approximations to control systems. Moreover, state 
observers are often used in steady state conditions where 
any disturbance is dampened or filtered by the system 
after transient conditions have been addressed. 
Measurement disturbance can be reduced and accounted 
for with the use of a reliable sensor with proven 
durability and low error range. With these system 
improvements the engine speed approximation can be 
ensured to stay within the desired percentage of the idle 
RPM that is considered the industry range for allowable 
speed deviations of ISC. Our work shows how state and 
measurement disturbance affects the error dynamics of a 
LTI state observer. Through simulation, we identified 
the effect of these disturbances and provided a thorough 
analysis of our findings. 
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Abstract subjected to a disturbance”. In power system environment 



In a deregulated power system structure, power 
producers and customers share a common 
transmission network for wheeling power from the 
point of generation to the point of consumption. All 
parties in this open access environment may try to 
produce the energy from the cheaper source for 
greater profit margin, which may lead to overloading 
and congestion of certain corridors of the transmission 
network. This may result in violation of line flow, 
voltage and stability limits. In power system 
environment voltage stability plays major role, it is 
integral part of the power system stability. In general 
Voltage stability problems occur more frequently in a 
heavily loaded system. Lack of voltage stability may 
develop severe problems in the power system, which 
might lead to system black out in the form of sudden 
voltage dips at the load buses. The cause of voltage 
instability could be numerous. The control of voltage 
levels can be accomplished using FACTS Controllers 
like Static Var Compensator (SVC), Thyristor-controlled 
phase angle regulator (TCPAR). This paper focuses on 
the enhancement Voltage Stability using FACTS devices. 
The optimal location of FACTS devices were determined 
based on Sensitivity methods. The effectiveness of 
proposed method is demonstrated on modified IEEE-9 
bus system using power world simulator 8.0. 

Keywords: Deregulation, Voltage Stability , TCPAR, 
SVC, Power world Simulator. 

1. Introduction 

In deregulation it is expected that the product will be 
exchanged at several points along the generation, 
transmission and distribution systems. Hence it is the 
responsibility of the power suppliers to supply safe and 
economical electric power to customers with the existing 
transmission line efficiently. 

“Voltage stability is the ability of a power system to 
maintain steady acceptable voltages at all buses in the 
system under normal operating conditions and after being 



voltage stability plays major role, it is integral part of the 
power system stability. In general Voltage stability 
problems occur more frequently in a heavily loaded 
system. 

2. Causes of Voltage Instability 

The main reason for voltage instability is the lack of 
sufficient reactive power in a system. Generator reactive 
power limits and reactive power requirements in 
transmission lines are the main causes of insufficient 
reactive power. 

2.1 Reactive Power Limits of Generators 

Synchronous generators are the main devices for voltage 
control and reactive power control in power systems. In 
voltage stability analysis active and reactive power 
capabilities of generators play an important role. The 
active power limits are due to the design of the turbine 
and the boiler. Therefore, active power limits are constant. 
Reactive power limits of generators are more 
complicated than active power limits. There are three 
different causes of reactive power limits that are; stator 
current, over-excitation current and under-excitation 
limits. The generator field current is limited by over- 
excitation limiter in order to avoid damage in field 
winding. In fact, reactive power limits are voltage 
dependent. However, in load flow programs they are 
taken to be constant in order to simplify analysis. 

2.2 Reactive Power Requirements in 
Transmission Lines 

Transfer of active and reactive power is provided by 
transmission lines. Since transmission lines are generally 
long, transfer of reactive power over these lines is very 
difficult due to significant amount of reactive power 
requirement. 

The characteristics of voltage stability and the effect of 
transmission lines are illustrated by 2-bus simple test 
system that is shown in Figure 1. 
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Active and reactive power delivered to load bus can 
be written as 

V V 

P = — — - Sind (1) 

X 

K V 2 „ _ V 2 

= — — - CosS (2) 

X X 



Eliminating S and solving Equation 1 and 2 for V- t yields 
Equation 



v 2 =A^--qx±x x 



4X 2 



P 2 -Q -fr (3) 

A 



In order to attain a real solution for V 2 Equation 3 should 
be satisfied. 



V 2 ~ V 2 

- l t - P 2 - Q — > 0 
4X 2 X 



(4) 



substitute the short circuit power at the receiving end, 

\ 2 



o = Yl_ (S ' 

sc x 



v 



J sc 

2 



> P 2 + Q S. 



SC 



(5) 



J 



The maximum possible transfer of active power is 



J sc 



(for Q — 0 ) and the maximum possible transfer of 



reactive power is 



J sc 



(for P= 0 ). Since 



J sc 



> 



J 



J sc 



V 



,it can be accomplished that transfer 



of reactive power is more complicated than transfer of 
active power and as it is observed from Equations 1 and 2 
transfer of power is inversely proportional to line 
reactance. 

Overhead lines either absorb or supply reactive power, 
depending on the load current. They generate reactive 
power under light load since their production due to the 
line shunt capacitance exceeds the reactive power 
requirement in the transmission line due to the line 
impedance. Under heavy load, they absorb more reactive 
power than they produce. 

Underground cables produce reactive power since the 
power requirements never exceed the production due to 



their high shunt capacitance under all operating 
conditions. Since they produce reactive power, a definite 
amount of reactive power demand of loads is met by the 
power produced by underground cables. Thus, the 
possibility of seeing voltage stability problem in buses 
decreases. 

3. Voltage Collapse 

Voltage collapse is the process by which the sequence of 
events accompanying voltage instability leads to a low 
unacceptable voltage profile in a significant part of 
system. 

When a power system is subjected to a sudden increase 
of reactive power demand, the required demand is met by 
the reactive power reserves supplied from generators and 
compensation devices. Most of the time, this can be 
achieved since there are sufficient reserves. Sometimes, it 
is not possible to meet this rapid increase in demand due 
to combination of events and system conditions. Thus, 
voltage collapse and a major breakdown of part or all of 
the system may occur. 

Typical scenario of a voltage collapse can be summarized 
as follows: 

• Some of the large generating units near load centers 
can be out of service due to abnormal operating 
conditions. This reduces reactive power supplied and 
some transmission lines are heavily loaded so as to 
carry sufficient reactive power to problematic parts 
of the grid. 

• Loss of a heavily loaded transmission line due to a 
fault causes additional loading on the remaining 
adjacent lines. This increases the reactive power 
requirement in lines since reactive power 
requirement increases rapidly for loads above surge 
impedance loading. Thus, reactive power demand in 
the system augments. 

• The load voltages decrease because of extra reactive 
power demand. This causes a decrease in load 
demand. The voltage control system of generators 
restores terminal voltage by increasing field 
excitation. The additional reactive power flow 
through transformers and transmission lines causes 
increased voltage drop across each of these elements. 

• As a result of tap changing operations, the reactive 
power output of generators increase. When the 
generator hits the reactive power capability limit, the 
terminal voltage decreases. Its sharing of reactive 
power is transferred to other generators, leading to 
overloading of more generators. The process 
eventually leads to voltage instability and voltage 
collapse. 

Therefore, it is very important to analyze the underlying 
principle of voltage stability in every power system. 

In assessing voltage security of power systems, we deal 
with P-V and Q-V curves. The methods of P-V and Q-V 
curves are widely accepted tools; they can provide real 
and reactive power margins, respectively. Variations in P 
and Q affect the voltages at the load buses. P-V and Q-V 
curves are generated by series of power- flow solutions. 
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3.1 Importance of PV Curve 

The process of voltage stability analysis involves the 
transfer of P from one region of a system to another, and 
monitoring the effects to the system voltages V. This type 
of analysis is commonly referred to as a P-V study. 

The P-V curves, real power- voltage curve, are used to 
determine the MW distance from the operating point to 
the critical voltage. A typical P-V curve is shown in 
Figure 2. 




Consider a single, constant power load connected through 
a transmission line to an infinite-bus. Let us consider the 
solution to the power-flow equations, where P, the real 
power of the load, is taken as a parameter that is slowly 
varied, and V is the voltage of the load bus. The three 
regions shown in Figure 2 are related to the parameter P. 
In the first region, the power flow has two distinct 
solutions for each choice of P; one is the desired stable 
voltage and the other is the unstable voltage. As P is 
increased, the system enters the second region, where the 
two solutions intersect to form one solution for P, which 
is the maximum. If P is further increased, the power- flow 
equations fail to have a solution. This process can be 
viewed as a bifurcation of the power-flow problem. The 
method of maximum power transfer determines critical 
limits on the load bus voltages, above which the system 
maintains steady-state operation. 

The P-V curve is drawn for the load bus and the 
maximum transmissible power is calculated. Each value 
of the transmissible power corresponds to a value of the 
voltage at the bus until V = Vcrit after which further 
increase in power results in deterioration of bus voltage. 
The top portion of the curve is acceptable operation 
whereas the bottom half is considered to be the 
worsening operation. The risk of voltage collapse is 
much lower if the bus voltage is further away, by an 
upper value, from the critical voltage corresponding to 
Pmax. Hence, the P-V curve can be used to determine the 
system’s critical operating voltage and collapse margin. 

3.2 Importance of QV-Curve 

The Q-V curves, reactive power- voltage curve, are used 
to determine the Mvar distance from the operating point 
to the critical voltage. A typical Q-V curve is shown in 
Figure 5.3. It shows the sensitivity and variation of bus 
voltages with respect to reactive power injections or 



absorptions. As a traditional solution in system planning 
and operation, the voltage level is used as an index of 
system voltage instability. If it exceeds the limit, reactive 
support is installed to improve voltage profiles. With 
such an action, voltage level can be maintained within 
acceptable limits under a wide range of MW loadings. In 
reality, voltage level may never decline below that limit 
as the system approaches its steady state stability limits. 
Consequently, voltage levels should not be used as a 
voltage collapse warning index. 




In Figure 3, the Q axis shows the reactive power that 
needs to be added or removed from the bus to maintain a 
given voltage at a given load. The reactive power margin 
is the Mvar distance from the operating point to the 
bottom of the curve. The curve can be used as an index 
for voltage instability. Near the nose of a Q-V curve, 
sensitivities get very large and then reverse sign. Also, it 
can be seen that the curve shows two possible values of 
voltage for the same value of power. 

The power system operated at lower voltage value would 
require very high current to produce the power. 

That is why the bottom portion of the curve is classified 
as an unstable region; the system cannot be operated, in 
steady state, in this region. The top portion of the curve 
represents the stability region while the bottom portion 
from the stability limit indicates the unstable operating 
region. It is preferred to keep the operating point far from 
the stability limit. 

The control of voltage levels is accomplished by 
controlling the production, absorption, and flow of 
reactive power at all levels in the system. The devices 
used for this purpose may be classified as follows 

• Sources or sinks of reactive power, such as shunt 
capacitors, shunt reactors, synchronous condensers, 
and static var compensators (SVCs) 

• Line reactance compensators, such as series 
capacitors. 

• Regulating transformers, such as tap-changing 
transformers and boosters. 

A static var system is ideally suited for applications 
requiring direct and rapid control of voltage. 

4. FACTS Controllers 

4.1 Static VAR Compensator (SVC) 

According to IEEE-CIGRE co-definition, a static var 
compensator is a static var generator whose output is 
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varied so as to maintain or control specific parameters 
(e.g. voltage or reactive power of bus) of the electric 
power system. 

SVC is a first generation FACTS controller that is 
already in operation at various places in the world. Static 
VAR Compensator (SVC) is a shunt connected FACTS 
controller whose main functionality is to regulate the 
voltage at a given bus by controlling its equivalent 
reactance as shown in Fig. 4. 




The reactive power injected into the bus due to SVC can 
be expressed as 

Q SVC = B SVC ^ ( 6 ) 



limit, the SVC becomes a shunt capacitor. At the 
inductive limit, the SVC becomes a shunt reactor. 



4.2 Thyristor Controlled Phase Angle 
Regulator 



TCPAR is identical with a phase shifting transformer 
with a thyristor type tap changer and plays an important 
role in increasing loadability of the existing system and 
controlling the congestion in the network. FACTS device 
like TCPAR can be used to regulate the power flow in 
the tie-lines of interconnected power system. 

This is also known as Static Phase Shifter (SPS) and 
phase shift with respect to the bus voltage is achieved by 
adding or subtracting a variable voltage component in 
quadrature with the bus voltage. The quadrature voltage 
is injected in series with the transmission line by a 
boosting transformer. The basic arrangement is shown in 
figure 4. 



\ m A\AJ 



TfVW 

y* = vi. 






Entiling Transformer 



Converter 



Boosting Transformer 

V* 

F s 



Where V t is the voltage magnitude of the bus at which 
the SVC is connected. 

This thesis provides information on the SVC total 
susceptance model required to achieve a given level of 
compensation. 



Figure4 Shows Static Phase Shifter 

5. Optimal Location Based on Sensitivity 
Approach for TCSC & SVC 




Figure .5 Voltage-Current Characteristics of SVC 



5.1 Optimal Location of Static Var Compensator 

For computing the loss sensitivity index with respect to 
SVC an exact loss formula has been used, which 
expresses Ql as, 



C 



dQt 

dQt 



2 Z"=1 (“ijQi+fijPj) i=1 > n ( § ) 



where a and B. are the loss coefficients defined as, 

ij v 



a.= 1 cos (8 i S.) 

J vy j 


(9) 


Pu= X ‘ j sm(S, 5.) 
1 vv 


(10) 



Figure .5 shows the steady-state and dynamic voltage- 
current characteristics of the SVC. In the active control 
range, current/susceptance and reactive power is varied to 
regulate voltage according to a slope (droop) 
characteristic. The slope value depends on the desired 
voltage regulation, the desired sharing of reactive power 
production between various sources, and other needs of 
the system. The slope is typically 1-5%. At the capacitive 



P + jQ is the complex injected power at bus i and 
V is the imaginary part of the ij th element of [Zbus]. 

The SVC should be placed at a bus i having most 
negative sensitivity index (C t ). The placement of SVC 
has been considered at load buses only. 
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5.2 Optimal Location of TCPAR 

Here, a method based on the sensitivity of the total 
system reactive power loss (Q L ) with respect to the 
control variables of the FACTS devices 

1. Loss sensitivity with respect to control 
parameter 0/ 7 - of TCPAR 
placed between buses i and j, 




dQi 

dO u 



2. Phase shift (ciy) for a TCPAR placed between buses 
i and j. In that line 0y is the net phase shift in the 
line and includes the effect of the TCPAR. The loss 
sensitivities with respect to Xy and 0y can be 
computed as: 

(“> 



5.3 Selection of Optimal Placement of FACTS 
Devices 

Using the loss sensitivities as computed in the previous 
section, the criteria for deciding device location might be 
stated as follows: 

l.SVC should be placed at a bus i having most 
negative sensitivity index 

2. TCPAR must be placed in the line having the highest 
absolute value of loss sensitivity index ( by ) 



6. Contingency 

The operation of a power grid requires unplanned 
component outages could lead to dramatic and costly 
blackouts. Planning and operating criteria are designed in 
such a way that “the interconnected power system shall 
be operated at all times as a result of any single 
contingency or multiple contingencies. 



6.1 Crticial Contingency Screening 

A suitable location of FACTS devices will be that which 
would give best performance in most of the operating 
conditions. To find out the severe line outage 
contingencies, a real power flow performance index (PI) 
has been used as defined in the following: 

f o \ 2n 



n = Y N ^ 

/ ^ w=1 In 



1 LM 
MAX 



P 

V 1 LM 



( 12 ) 



Where P LM is the real power flow 

P^ X is the rated capacity of line-m, 
n is the exponent and 

W m a real nonnegative weighting coefficient 
which may be used to reflect the importance of lines. 



PI will be small when all the lines are within their limits 
and reach a high value when there are overloads. Thus, it 
provides a good measure of severity of the line overloads 
for given state of the power system. By most of the 
operational standards, the system with one huge violation 
is much more severe than that with many small violations. 



In this study, the value of exponent has been taken as 2 
and Wi =1. 



7. Results & Discussions 

The Proposed methods for optimal placement of TCPAR 
and SVC for voltage stability improvement are tested on 
IEEE 9 bus system. The Computation of power flow is 
carried out by using Power world simulator to calculate 
the bus voltages. The Critical contingency screening is 
performed using Performance index using MATLAB 
programming. 

7.1 System Model 

This system consists of 9 buses, 9 line sections, 3 
generator buses and 3 load buses. 




Figure .8 Single line diagram of IEEE 9 bus System 

7.2 Sensitivity Index 



Table 1. Var loss sensitivity index of IEEE - 9 bus system 



Line 
No. (K) 


From- To Bus 


Loss Sensitivity 






TCPAR 

(bn) 


SVC (Cy) 


1 


1-4 


1.441 


-3.4022 


2 


8-2 


-3.26 


-3.438 


3 


3-6 


1.696 


-3.263 


4 


4-5 


0.588 


4.1572 


5 


9-4 


-0.752 


0.4279 


6 


5-6 


-1.19 


0.9455 


7 


6-7 


0.4486 


-2.4924 


8 


7-8 


-1.513 


-1.6932 


9 


8-9 


1.679 


-3.9318 



For this system, from Table 1, 

1. Optimal location of TCPAR for an IEEE 9 bus 
system based on loss sensitivity analysis the line 
having the highest absolute value of loss sensitivity 
index by index is chosen for placement of the 
TCPAR device. Hence bus 3 is selected. 

2. Optimal location of SVC for an IEEE 9 bus system 
based on loss sensitivity analysis is the lines having 
the highest negative sensitivity index is chosen for 
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placement of the SVC device. Hence bus 9, a load 
bus is selected. 

The optimal locations considered are line3 (3-6) for 
TCPAR with a phase shift of 2.9 degrees and unity 
tap ratio and bus 9 for SVC with a nominal value of 
50 Mvar. 

Modified IEEE 9bus system incorparating 
TCAPAR,SVC shown in 
Figure. 9 




Figure. 9Single line diagram of IEEE 9 bus System with FACTS devises 




Bus Number 



■ Without FACTS ■ With TCPAR a With SVC 



Figure. 10 Comparison of voltages for 9 bus system at base case 

Figure 10 shows the comparision of voltages for an IEEE 
9 bus sytem without and with FACTS devices under 
normal conditions. It is observed that by placing SVC at 
bus 9 the voltage is more improved than by placing 
TCPAR line3 (3-6) 

The critical contingency (line outage) ranking using the 
PI values, are computed for each single line outage case. 
The ten critical lines at the base case loading, in the order 
of relative severity are given in the Table. 2 By observing 
the Table .2, line 9(8-9) is the most critical line for 
outage as contingency which is having the highest PI 
index. 

The voltage profile at each bus by placing TCSC and 
SVC for single outage of line at line 9(8-9) is shown in 
Fig. 10. It is observed that by placing SVC at bus 9 the 
voltage is more improved than by placing TCPAR line3 
(3-6)even at single outage of line. 



Table 2. Single Outage Contingency Ranking for 
IEEE 9 Bus System 



S.No 


Line No 


End Buses 


PI Index 


1 


Intact Case 


- 


0.504 


2 


1 


1-4 


0.86 


3 


2 


8-2 


0.87 


4 


3 


3-6 


0.97 


5 


4 


4-5 


0.81 


6 


5 


9-4 


0.83 


7 


6 


5-6 


1.41 


8 


7 


6-7 


0.69 


9 


8 


7-8 


0.844 




1 2 3 4 5 6 7 8 9 

Bus Number 

■ Without FACTS ■ With TCPAR ■ With SVC 



Fig ure.10 Comparison of voltages for 9 bus system at single outage 
(8-9) 

Table 3. Double Line Outage Contingency Ranking for IEEE 9 Bus 
Jllfstem 



S.No 


Line No 


End Buses 


PI Index 


1 


Intact Case 


- 


0.504 


2 


9,6 


8-9, 5-6 


1.24 


3 


9,3 


8-9, 3-6 


0.974 


4 


9,2 


8-9, 8-2 


0.884 


5 


6,3 


5-6, 3-6 


0.67 


6 


6,2 


5-6, 8-2 


0.942 


7 


3,2 


3-6, 8-2 


1.72 




1 2 3 4 5 6 7 8 9 

Bus Number 

■ Without FACTS ■ With TCPAR I With SVC 



Figured 1 Comparison of voltages for 9 bus system at double outage 
(3-6, 8-2) 
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From Table 3, the top five locations are selected and PI 
values are calculated for the different combinations of top 
five locations to find the severity of the system with 
double outages. Table 3 shows the double line outage 
contingency ranking for an IEEE 9 bus system. By 
observing the Table 3, line 3(3-6) and line 2(8-2) are the 
most critical lines for double outages as contingency 
having the highest PI index. 

The voltage profile at each bus by placing TCPAR and 
SVC for double outage of lines at line 3(3-6) and line 
2(8-2) is shown in Fig.ll. It is observed that by placing 
SVC at bus 9 the voltage is more improved than by 
placing TCPAR in line (3-6) even at double outage of 
line. 




Fig 12. shows the PV curve at bus 5 for an IEEE 9 bus 
system without and with FACTS Devices. From the 
Fig. 5. 9, it is observed that the voltage stability margin 
has been more improved by placing SVC at bus 9 than by 
placing TCPAR line3 (3-6). 




Figure. 13 QV Curve at bus 5 for 9 bus system 

Figure 13. shows the QV curve at bus 5 for an IEEE 9 
bus system without and with FACTS Devices. From the 
figure 13, it is observed that the voltage stability margin 
has been more improved by placing SVC at bus 9 than by 
placing TCPAR line3 (3-6). 



Lack of voltage stability may develop severe problems in 
the power system, which might lead to system black out 
in the form of sudden voltage dips at the load buses. The 
cause of voltage instability could be numerous. The 
control of voltage levels can be accomplished using 
FACTS Controllers like Static Var Compensator. 

It is observed that the SVC consistently gives the better 
performance than the TCPAR in all case studies 
conducted, also giving better voltage stability margin. 
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8 Conclusion 

In power system environment voltage stability plays 
major role, it is integral part of the power system stability. 
In general Voltage stability problems occur more 
frequently in a heavily loaded system. 
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Abstract 

This paper presents inrush current phenomenon in a 
single phase and 3 -Phase transformer. This is initially 
done in the absence of any inrush current suppressing 
device and then with Pulse Width Modulated (PWM) 
voltage source converter, Finally it is compared with 
other technique having DC voltage regulation using PI 
controller. In both the proposed techniques main 
transformer is connected in series to a PWM voltage 
source converter through a matching transformer to the 
source. As the PWM converter is connected in series to 
the source it acts as a resistor for the source current there 
by producing no inrush phenomenon. Inrush current 
suppression principles for both methods were discussed 
with the help of block diagrams. Variation of 
magnetizing current with varying control gains is 
discussed in both cases and the results were shown to 
confirm the validity. Both methods were developed & 
simulated in MATLAB/SIMULINK environment. 1 

Keywords : Inrush current, PWM converter, PI 

Controller, DC regulation. 

Nomenclature 

V re f Reference output signal for converter 
K Control gain 

I s Source current 

Vt r ab Main transformer line to line voltage 

V 0 Voltage across converter 

Vtr Voltage across transformer 

I sa ‘a’ phase source current 

V oa Voltage across PWM converter for ‘a’ phase 

1. Introduction 

Switching of the transformers will be done in many 
times during a day especially when the transformers are 
used in industrial facilities. It is well known that each 



1 This study has been implemented and both the methods 
were developed & simulated in MATLAB/SIMULINK 
environment 



energization of the transformer will generate a transient 
inrush current even if the transformer carry no load. 
These inrush currents cause mal-operation of the 
protective relays which needs to be avoided. If it is a 
single phase transformer, inrush phenomenon can be 
avoided easily by simply connecting the transformer to 
the source at closing angle of 90 or 270 degrees. But, in 
3 -phase circuits, the inrush phenomenon cannot be 
avoided by simply adjusting the closing angle since the 
phase angles of the source voltages differ from each 
other by 120 degrees. Inrush phenomenon can be 
avoided by connecting series resistors which needs to be 
shorted before loading the transformer on the secondary 
side. But this method may cause re-inrush phenomenon 
with increased losses [1]. A PWM Converter can be 
connected in parallel which acts as a inrush current 
sinker on the source side but requires a large rating for 
the PWM Converter [2]. 

In the present work, proposed method of Suppressing 
the magnetizing current with series connected voltage 
source PWM converter is verified and extension of the 
method with dc voltage regulation using PI controller 
was discussed. With series connected voltage-source 
PWM converter the required rating of the PWM 
converter is about 0.25% of the main transformer in 
single phase circuits and it is less than 0.1 % in 3 -phase 
circuits illustrating that it is useful for practical 
applications [3]. 

The rest of the paper is organized as follows; section 2 
explains the inrush current of transformers. Section 3 
presents suppression of inrush current using PWM 
converter. Section 4 introduces inrush current 
suppression with DC voltage regulation. Finally section 
5 is the conclusion. 
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2. INRUSH CURRENT OF 
TRANSFORMERS 



Figure 1 Constructed transformer model for 1 -phase. 

A 1 -phase 200V, 60Hz, 6KVA base transformer is 
connected to an AC voltage source consisting of 
resistive impedance through a breaker & current 
transducer as shown in Figure 1. Breaker implements a 
Circuit breaker opening depending upon the external 
control of switching times. 

In the present work a resistive snubber is considered. 
With the help of the current transducer, inrush current is 
measured and the results were shown for both closing 
angels of 0 degrees & 90 degrees. It was observed that 
high inrush current (5 3 A) occurred when the transformer 
is connected with a closing angle of zero degrees & the 
inrush current is having less magnitude (0.5 A) when the 
transformer is connected with a closing angle of 90 
degrees. Figure 2 and Figure 3 represents Inrush Current 
of single phase transformer with closing angle of zero 
and 90 degrees respectively. 




Time 

Figure 2. Inrush Current of transformer with closing angle of 0° 




Figure 3. Inrush Current of transformer with closing angle of 90° 

The same basic model is extended for 3 -Phase circuits 
and the source current waveforms on the primary side of 
the 3 -phase transformer are shown. As the voltages of 



all the 3 -phases differ with an angle of 120 degrees, with 
closing angle adjustment, we cannot avoid inrush 
current. It was observed that even with closing angle of 
90 degrees, there is a flow of inrush current (120 Amps) 
in all the phases unlike single phase transformer circuits, 
where it was observed as 0.5 Amp’ 




Figure 4. Three Phase Inrush Currents of transformer with closing 
angle of 90 degrees 

3. SUPRESSION OF INRUSH CURRENT 
USING PWM CONVERTER 

1. A series voltage source PWM converter is connected 
through a matching transformer between source and 
main transformer where actual inrush phenomenon 
occurs. 

2. In order to suppress the switching ripples caused by 
switching action, on ac side of the converter LC filter is 
connected. 

3. The Control gain of PWM converter will be 
multiplied by the primary side source current of the 
transformer in order to obtain the reference output signal 
for the converter which means voltage source PWM 
converter is acting as a resistor of K[Q] for the source 
current is I s . 

V ref = K.I S (1) 

4. In the present work, for generation of gate pulses to 
PWM converter, sine triangular intercept technique is 
used. The primary side source current after passing 
through the low pass filter is treated with the triangular 
wave depending upon the required variation of control 
gain of PWM Converter. 

Initially with the help of switch 1 main transformer is 
connected to the source at zero output voltage. Voltage 
source PWM converter will act as a resistor in order to 
suppress the inrush current. Finally after suppressing 
inrush current, switch 2 will be closed which results zero 
output voltage for PWM converter that is PWM 
converter is disconnected from the main circuit. 

Even at closing of Zero or 180 degree where high 
magnitudes of inrush current results, we have to design 
in a manner such that the source current does not exceed 
the open circuit current. With open circuit current & 
voltage of the transformer we can get resistance value 
that is required for the PWM converter. The control gain 
K is determined as 765Q considering modulation index 



Breaker 




Current 




Saturable 






Transducer 




Transformer 



Ac voltage Source 
with Impedance 
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of 0.8. Data in Table-2 represents the values of the 
circuit constants used in calculations & simulations. 



Turn Ratio 


1:1 


Open Circuit Current 


0.2A 


Open Circuit Power 


25W 


Short Circuit Power 


1 16W 


Source Voltage 


200V 


Source Frequency 


60Hz 


Control gain 


765 Q 


DC Capacitance 


500jaF 


Initial Vol. of DC capacitance 


270V 



TABLE I: DATA OF l-$ 200V, 50 HZ, 6KVA 
TRANSFORMER 

In 3 -Phase circuits, three individual circuits will be 
connected as shown in Figure 6. The same concept of 
inrush current suppression is used by placing 3 -phase 
PWM converter between the source voltage and the 3- 
phase transformer. Individual gains will be calculated 
and as the large control gain K of the voltage source 
PWM connecter is effective for reducing the settling 




Switch 2 



time of magnetic flux of core, large control gain of 
765Q is used in the present work to reduce the operating 
time. Longer operating times results larger amount of 
active power flow into the dc capacitor. 

Single phase and three phase circuits were observed for 
different variations of control gain of PWM converter. 
With control gain suddenly made Zero i.e. removing the 
pulses to the converter suddenly at a particular instant is 
leading to re- inrush current phenomenon, with re-inrush 
current magnitude of 2.45 Amps for a while then to 0.35 
amps at the primary or source side. Here in this case, the 
pulses are supplied till 0.15 sec and after 0.15 sec, 
converter was switched off suddenly. It was also 
observed that at 0.15 sec, even after removing the pulses 
a small voltage of 5V was observed across switch 2 
which may be because of breaker resistance. Figure 7 
Illustrate the simulated waveforms when the voltage 
source PWM converter is suddenly stopped in single 
phase circuits. In all the cases V tr represents the voltage 
across the transformer, I s represents the source current of 
the transformer and V 0 represents the voltage across 
PWM converter. 




Figure 5. Inrush Current suppressor with PWM converter in 1- phase circuits. 
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Delta 




Figure 6. Inrush Current suppressor with PWM converter in 3- phase circuits 






Time 

Figure 7. Re-inrush phenomenon when suddenly K made zero 

Figure 8 represents the waveforms when the control gain 
of the PWM converter after 0.15 sec is gradually 
decreased over a period of 40ms from its initial constant 
value using ramp signal prevents re-inrush phenomenon 
and source current is observed with a magnitude of 0.35 
Amps. It was also observed that with gradual decrement of 
control gain before stopping the voltage source PWM 
converter results zero output voltage across converter in 
addition to the suppression of re-inrush phenomenon. 




Figure 8. l-(j) output waveform’s with gradual decrement of control gain 

Inrush current suppressor is applied for 3 -Phase circuits, 
and the simulation waveforms are shown in Figure 9. 
with gradual decrement of control gain of voltage source 
PWM Converter, ‘a’ phase voltages and currents are 
represented in Figure 9. V tra b is the line to line voltage 
between a & b phases at the main transformer terminals. 
I sa is the ‘a’ phase primary side current and V oa is the 
voltage across the PWM Converter. We observed that the 
series connected PWM converter suppressed the inrush 
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current perfectly in both single phase & 3 -phase circuits .It 
was also observed that without series voltage source 
PMW converter active power is high and with PWM 
converter it was low. 




0 0.05 0.1 0.15 0.2 0.25 0.3 




Time 



Figure 9. Waveform’s with gradual decrement of control gain in 3-Phase 
Circuits. 



4. INRUSH CURRENT SUPPRESSION 
WITH DC VOLTAGE REGULATION 

For the same circuit another method is proposed to reduce 
the inrush current by using dc voltage regulation. Pulses 
are generated considering the voltage across the capacitor 
of the converter. By treating V dc with low pass filter and 
using PI controller the pulses are generated. Control gain 
of the converter is considered as 765 and the initial 
voltage of DC capacitor as 270V. Figure 10 represents 
the block diagram approach for pulse generation using V dc . 

In order to match the digital & analog signals an integrator 
used with these pulses. It was observed that with gradual 
decrement of control gain, about 40 ms after 0.15 sec, the 
magnitude of inrush current was 0.15 Amps as shown 
Figure 1 1 where as with the first method it was 0.45 Amps. 




Figure 10 .Block Diagram for pulse generation using V dc 




Figure 1 1 . Inrush current with gradual decrement of K 

When the pulses were suddenly stopped at 0.15 sec, a 
re-inrush phenomenon was observed with re-inrush 
current magnitude of 1 .05 Amp and magnetizing current 
magnitude of 0.15 Amps. Figure 12 represents the 
primary current when K is suddenly stopped. 
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Figure 12. Inrush current when K suddenly stopped 

By using the first method, re-inrush current magnitude 
is 2.45 Amps & inrush current magnitude 0.35 Amps, 
shows that proposed method is having best results than 
first method. 
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lime 

Figure 13. Output wave forms for 3 -phase circuits with gradual 
decrement of K. 

Figure 13 shows waveforms for 3 -phase circuits with 
DC regulation and gradual decrement of control gain of 
voltage source PWM converter. Phase voltages and 
currents are shown in Figure 13. V tr is the line to line 
voltage between a & b phases at the transformer 
terminals. I s is inrush current per phase and Vo is the 
voltage across the series connected PWM Converter. 
Without DC voltage regulation primary side source 
current was observed as 0.35 Amps and with dc voltage 
regulation it was observed as 0.15 Amps. Reactive and 
active power waves forms also observed for all the cases. 

5. CONCLUSIONS 

In this paper, the confirmation of the inrush 
phenomenon was done for a single and 3 -phase 
transformer models using MATALAB/SIMLINK 
environment. It was observed that by changing the 
closing angle of circuit breaker we can control the 
magnetizing current in single phase transformers and it 
is not possible for the 3-phase transformers. Validation 
of inrush current suppression with series connected 
PWM convertor was done with variations in control gain 
of convertor in both single and 3- phase circuits. It was 
observed that with gradual decrement in control gain 
before stopping the pulses to series voltage source PWM 
converter results zero output voltage across the 
converter in addition to the suppression of re-inrush 
phenomenon. 

A new technique was proposed to suppress the inrush 
current with DC voltage across the capacitor. With DC 
voltage regulation and variation in control gains there is 
a further decrement in the magnitude of inrush current 
which shows that it is more effective technique than 
series connected voltage source PWM convertor 
technique. The proposed method applied to 3 -phase 
circuits also and satisfactory results were obtained. 



Active and reactive power for all the cases were 

observed. It can be concluded that the proposed method 

is superior than the method using series connected 

voltage source PWM convertor. 
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ABSTRACT 

Although the modern torque converter (TC) has many 
functions, this paper focuses on its ability to transmit 
torque while allowing the engine and gearbox to slip 
or rotate at different speeds. While this slipping 
greatly improves drivetrain feel, it can also reduce 
driveline efficiency due to friction. To improve 
efficiency, modem torque converters include a bypass 
Torque Converter Clutch (TCC), which allows the 
slipping to be reduced or eliminated at various 
operating conditions. This paper will investigate how 
to apply Piecewise Linearization to the plant model, 
so it can be controlled by switched gain controllers. 
These controllers are tuned by both Lyapunov and 
Riccati methods, and their performance evaluated 
against the conventional PID controlled model in 
terms of settling time, response and performance 
robustness. 

Keywords: Lyapunov; Riccati; Control theory; 
Torque Convert Clutch; Piecewise Linear; 
Measurement disturbance; Estimation methods. 

Nomenclature: TC - Torque Converter; TCC - 

Torque Converter Clutch; AMT - Automated Manual 
Transmission; DCT - Dual Clutch Transmission; PID 

- Proportional Integrated Differential; PWL - 
Piecewise Linearization; State matrixes A, B, C, D; X 

- system states; Y - system outputs; R, Q and y - 
tuning parameters for Riccati method 

1. INTRODUCTION 

Automatic transmissions have been the choice of 
many drivers for their smooth operation since their 
introduction. Much of their smooth gear changes and 
isolation from engine vibrations is accomplished by 
fluid coupling in the torque converter between the 
engine and the gearbox. With increasing demand for 
vehicle gas mileage and efficiently, much attention 
has been directed to improving driveline efficiencies. 
One area to improve driveline efficiency is in 



automatic transmissions. To address these new 
demands on automotive transmissions, modem 
vehicles have begun to adopt new types of 
technology, such as automated manual transmission 
(AMT), a subclass of AMTs the dual clutch operated 
transmission (DCT), and the constant velocity 
transmission (CVT) into the marketplace alongside 
traditional hydro-mechanical fixed-ratio gear 
gearboxes with TC. Although transmission 
efficiencies in AMT have been greatly increased by 
replacing TCs with friction clutches, their smoothness 
has gone down significantly because of the 
quantitative nature of friction clutch engagement. 
Much work was been done to control TC engagement 
by improving the driveline feel, but AMT are still no 
match for traditional TC equipped vehicles. 
[1,2, 3, 4, 5] The objective of the powertrain designers 
is to control the drive clutch engagement to mimic the 
smooth and overdamped response of TC equipped 
vehicle with the underdamped nonlinear system 
behavior. Fuzzy Logic, optimal control and Kalman 
Filters control schemes have been used in an attempt 
to replace look up tables to achieve the desired 
smooth underdamped response however this paper 
will investigate using switched gain linear controllers. 
[6,7,8] CVTs, on the other hand, are known for their 
ultra- smooth underdamped operation, but their lack 
of torque capacity, “mush box” reputation and high 
cost have limited their penetration into the market 
place. [9] 

To keep costs low and to maintain the current 
automatic transmission shift quality, most 
manufacturers have chosen to improve the efficiency 
of the more common torque converter/fixed ratio 
transmission found in most new vehicles sold in the 
United States by improving the TCC control. This 
paper discusses the advantages of improved dynamic 
modeling and control of these transmissions. 
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In fixed-ratio automatic transmissions, one of the 
most important components is the torque converter 
and its fluid coupling between the engine output and 
the gear box input. This fluid coupling is necessary 
to allow the engine and transmission to turn at 
different speeds when the transmission is changing 
gears or when the vehicle is stopped. A secondary 
benefit of this slip, the speed differential between the 
clutch’s input and output, is that it multiplies the 
engine torque as the slip speed increases. A third and 
less studied benefit of this slip is that it provides 
engine isolation from the drive train by absorbing the 
engine’s rotational accelerations, or roughness, at low 
engine speeds. 

Like all systems, there are tradeoffs; torque 
converters reduce driveline efficiency by diverting 
input power into heating the transmission fluid 
whenever the torque converter is slipping. TCCs 
reduce the slipping by bypassing the TC with a 
mechanical clutch as shown in figure 2. Although the 
TCCs improve efficiency, they do it at a cost: 
increased driveline harshness. [10,11] This 
harshness comes in two primary ways; the jerk 1) 
from locking up the clutch and 2) by transmitting 
underdamped engine vibrations. This paper focuses 
on the latter and how vibrations at lower engine 
RPMs are transmitted to the rest of the driveline. 
Although locking up the clutch improves efficiencies 
by eliminating torque converter slip and allowing the 
engine to run at a lower speed, this reduced internal 
combustion engine (ICE) often creates rotation speed 
fluctuations. These speed fluctuations are created 
when the individual cylinders ignite during their 
power stroke creating torque impulses. On a 4 
cylinder four stroke ICE, there are two power strokes 
per revolution. With a typical lockup speed around 
1500 RPM, this is a noticeable 50 Hz/second torque 
impulse causing fluctuations in engine speed. These 
speed fluctuations are dampened by the TC and to a 
lesser extent other driveline components. However, 
when the TC is bypassed, much of it is directly 
transmitted to the drivetrain, which increases vehicle 
vibration and noise. [10,12] 

Clutch slipping control is a technique in which a 
partially engaged TCC is allowed to slip a limited 
amount thus preserving the damping between the ICE 



and gear box while retaining most of the benefits of a 
locked up TCC. By keeping the slip to one or two 
hundred RPM, the ICE’s speed variations can be 
absorbed by the slipping TCC, allowing it to be 
applied in a wider range of conditions while still 
meeting comfort objectives. 

While the mechanics of how clutches slip and how 
they are controlled in this region of operation are well 
documented [13,14,15,16,17], this paper discusses 
how Pricewise Linear (PWL) plant models are used 
to design switched gain controllers. These controllers 
are then optimized to create faster and more robust 
systems than the current nonlinear models controlled 
with PID controllers. 

This paper discusses piecewise linearization in 
section 2 and applies it to the nonlinear SimuLink 
model in section 3. Section 4 presents the baseline 
results compared with a new State controlled 
switched gain controller, shown in figure 1, for its 
settling time and robustness to model speed 
measurement noise. The new controller was 
optimized using two different methods, Lyapunov 
and Riccati, and their results compared to baseline 
model results. Section 5 presents the analysis of the 
findings and conclusions. 
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lock-up 

torque converter 



fluid pressure 
applies clutch 



agcoauto.com 



Figure 2 - Torque Converter Clutch Diagram 



2. PIECEWISE LINEARIZATION 
(PWL) OF MIMO MODELS 

One of the most difficult aspects of modeling 
complex, multi input multi output (MIMO) systems, 
such as an automobile’s drive train, is classifying the 
relevant system states and the individual subsystem 
transfer functions. This paper addresses these 
challenges by creating a PWL model of the system 
and implements a controller for each segment’s range 
of operation. Using the Lijun/Nianjiong/Daojun [10] 
description of the slipping control area, we can make 
assumptions that throttle input and vehicle speed 
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remain relatively constant in the slipping control area 
while clutch control input nonlinearly affects the 
amount of slip. [11] For simplicity, this paper does 
not address other slowly changing dynamic factors, 
such as fluid heat, clutch wear and vehicle load, 
which are dealt with in other papers. [18,19,20] The 
system in question can be described as follows: 

x(t) = f(x,u) (1) 

y(t) = h(x, u), 

where the input is u (clutch control input), the states 
are TCC input acceleration, or engine acceleration, 
and TCC output acceleration, or transmission input 
acceleration, and are defined as x and the output of 
the system is y (TCC input speed and TCC output 
speed). 

Then states x and output y are approximated by PWL 
estimates and represented by x and y, respectively. 
These estimates are calculated by choosing an 
operating point of the function, finding its 1 st 
derivative and using these to generate an estimate 
using the following form: 
x — f'(x ) * x + b. 




X 

Figure 3 - PWL Estimation of f(x) with 1 Segment 



When this is expanded into segments, it is 
represented by: 
x = fi(x) * x + b t 

where it is indexed by i = 1, 2, 3, ... , q for each 
segment and q is the number of segments. This 
method models the nonlinear function 
f(x 1 (t),x 2 (t), ...,x n (t)) as a collection of linear 
functions in the form: 
f'(x 1 (t),x 2 (t),...,x n (t)) = m ll x 1 {t) + 

b l m l2 x 2 (t) + ... + m in x n {t) + b t . 
The resulting PWL state space equations are: 

x(t) = A{x(t ) + B t u(t ) and (2) 

y{t) = C t x(t) + DiU(t), (3) 

with Ai, Bi, Q and D[ representing the system, input, 
output and feed forward matrices. 

An ideal example is the nonlinear electrical energy 
required to move the hydraulic solenoid actuating the 
clutch. Knowing that W = i 2 * R, where W the 
electrical power in watts, i is current and R is 
electrical resistance of the solenoid, it can be seen 
that the clutch force is nonlinearly proportional to the 
current. In this example, the nonlinear part of the 



power equation is a function of /(£) = i 2 . By 
inspection of /(£), it is apparent that if /(£) is 
linearized around i = 1, then a good estimate is 
/(£) = i. However, if it is linearized around i = 10, 
then the reader quickly sees /(£) = £ * 10. 
Depending upon many factors, including application, 
observability of the function, noise, and robustness of 
the control, this error can quickly destabilize the 
controller. To reduce the error, the function can be 
broken into segments with each segment having its 
own linear model. Even with little practical 
experience, it is clear that as the number of segments, 
q, increases, the error of / (£) becomes smaller. 
Following the processes described by Lowe (5), 
fi(i) = m{x{t ) + b t with i = (4) 

will be linearized from 0 to 100. Some rules that 
must be observed during PWL estimations are: 

• The function must be differentiable. 

• The center point of each segment must be found. 

• Discontinuities should be avoided. 

• A larger value of q equates to smaller estimate 
errors. 

Example 1 with q=l 

With only 1 segment, the center point is at t = 50. 
Differentiating /(£)' = 2 * i and substituting 
i = 50 gives /(£)' = 100. Then calculate the value 
of f at the center point: /( 50) = 50 2 . The error we 
want to minimize is defined as error = | /(£) = 
/(£) |. From (4) /(£) = /(£)' * x(i) + b and 

when solved for i=50, b = - 2500. The linear estimate 
with q=l of /(£) is as follows: 

/(0 = 100 *x(Q - 2500. (5) 

The results of the nonlinear response vs. the PWL 
model with one segment is shown in Figure 3. 




X 

Figure 4 - PWL Estimation of f(x) with 2 Segments 



Example 2 with q=2 

Similarly, if the t segment is divided into two 
segments, the center points are 25 and 75 yielding 
A(0 = 50*x(Q - 625 for 0<= i <= 50 

and 

/ 2 (i) = 150*x(Q - 5625 for 50< i <=100 
and is plotted as follows: 

It can be shown that as q approaches infinity, the 
linearization error approaches 0 for all values of i. In 
the next section, apply this to the Torque Converter 




27 




ACSE Journal, ISSN: Print 1687-481 1 , Online 1687-482X, CD-ROM 1687-4838, ICGST LLC, Delaware, USA, December, 2014 



Clutch model. The results of the nonlinear response 
vs. the PWL model with two segment is shown in 
Figure 4. 

3. PIECEWISE LINEAIZATION OF 
CLUTCH MODELS 

This paper’s nonlinear TCC is modeled in 
MatLAB/Simulink in Figure 5. The Inputs are engine 
torque request and clutch engagement control, while 
the outputs are engine output speed and transmission 
input speed. From these outputs, clutch delta, vehicle 
speed and vehicle acceleration can be easily 
calculated. 




Figure 5 - Simulink Model of Nonlinear Drive train 

For simplicity, this paper holds the throttle input 
constant at 4.5%. Because the nonlinear transfer 
function of the clutch control, system input, to the 
slip speed, system output, is not fully known, the 
TCC input and outputs at various steady state where 
recorded to create the data map in Table 1 . This table 
plots the normalized Clutch Input (Column 1), 
Throttle Input (Column 2), TCC Input speed (Column 
3), TCC Output speed (Column 4), and the assigned 
Segment number for the data range TCC Input speed 
(Column 5). These 35 data points were broken into 
operational segments and numbered 1 through 7. It 
should be noted that the difference between the 
engine speed and transmission input speed is the slip 
speed and is controlled by the clutch, creating a single 
input two output system. 

The typical engine speed, TCC input, and 
transmission input speed, TCC output, are plotted 
against the normalized clutch control, the system 
input, in Figure 6. 

For computation efficiencies and to exploit the 
common linear techniques, many methods have been 
used to estimate the response of these systems 
linearly [21,22]. This paper approximates the 
nonlinear function by dividing the function into linear 
segments using the Lowe/Zohdy method [23]. The 
superiority of this method over others is that any 
linear control design method may be used as long as 
it is continuous, thus eliminating the necessity of 
finding derivatives. By setting the output of the 

system, y(t) , equal to X(t), C = [J and 



D = [jjj and x(t) = /(u(t)) = m * u + b results 

in the new state equations: 

x(t) = A t * (nii * u + b t ) + B^u 

and y(t) = ra* * w + (6) 

2000 
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1700 
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^^TCC Input (Engine Speed) 
^^TCC Output (Gearbox Speed) 

Figure 6 - Normalized Clutch Input vs. TCC Speeds (RPM) 

where A t and B t are the PWL segmented state and 
input matrices segmented by i = 1,2,3, ... ,q, where 
q is the number of segments, and ra* * u + b t is the 
segmented linear approximation of x as a function of 
u. Where f (x ± (t) , x 2 (t) , ... ,x n (t )) is a collection of 
linear functions and takes the form: 
f(x 1 (t),x 2 (t),...,x n (tj) = m il x 1 (t) + 
bi m i2 x 2 (t) + ... + ra in x n {t) + b t where f, m* and 
bi are independent for both the engine function and 
the transmission function yielding 

f engined) ~ *n engine i * ^ b erL gi ne . 

and 

ftrans GO — ^4 trans j ^ T bt ra ns i 

i = 1,2, ...,q. 

Observing the segmenting guidelines described in the 
Lowe/Zohdy paper, the nonlinear transfer function 
was segmented into 7 pieces indexed by q and 
linearized using the least squares error, LSE, method. 
When the original engine and transmission speed 
derivatives, x engine and x trans , are overlade with 
their PWL estimations it is shown that the RPM error 
is within an acceptable few percent over most of the 
operating range of clutch control inputs in Figure 7. 
By linearizing the speeds in figure 1 as x = m x * u + 
b x and substituting into equation (2), a new equation 
for their acceleration is: 

x = A t * ( m Xi * u + b x . ) + Bi*u. (7) 

Setting this equal to the PWL estimations of the 
acculturations in figure 2, x = m x * u + b x , the state 
space matrix A t and B t can be solved algebraically 
as: 
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At = ^ and B t = m k - (L) * m x . (8) 

Using these State matrixes, a PWL model of the 
nonlinear drive trains was created. The results of 
both the nonlinear system and this model controlled 
by the original PID controller are shown in figure 8. 




TCC Input • TCC Output 



By setting (8) and (9) equal, we have derived a new 
closed-loop model system 

x(t) = A m x with A m - A- Bk . (10) 

In practice, it is often advantageous to calculate the 
open loop system poles, the eigenvalues of A, to 
determine its stability, and then choose the closed 
loop system poles, the eigenvalues of A m ,to be 2 to 4 
times faster, making sure they are all <= 0 to ensure 
closed loop stability. 

The system’s stability is ensured by the following 
Lyapunov equation: 

PA m + A m T P = -Q 

Given any Q, a symmetric positive-definite matrix, 
there exists a symmetric positive-definite matrix P. 
The solution for finding P can be found in the 
literature [27,28]. 

The simulation results of this method to a 100 delta 
slip speed request are shown in figure 9. 



#• PWL TCC Input •• PWL TCC Output 



Figure 7 - Normalized Clutch Input vs. 
Nonlinear and PWL TCC Speeds (RPM) 
PID Controlled Nonlinear and PWL Model 




Figure 8 - Nonlinear and PWL Delta Speed Request Response 
w/PID Controller 

LYAPUNOV OPTIMIZED CONTROLLER 

Lyapunov optimized control schemes have been the 
primary technique used for controlling nonlinear 
systems since its adoption to computers over 30 years 
ago. [24,25,26] The basic principle behind 
Lyapunov’s closed loop control law is that if the 
system can be defined in such a way that its energy is 
decreasing, a controller can be designed to make it 
stable at all times. This hinges on defining the energy 
as a positive and forcing its derivative to be negative 
definite or negative semi-defmite thus ensuring the 
values are positive, decreasing and limited, creating a 
asymptotically stable response to any input. 
Expanding on (1) above, we are replacing the input 
with u = u(x), giving the following: 

x(t) = f(x,u(x)). (6) 

By restricting u(x) to scalar gain k of the state, the 
input becomes 

u = -kx (7) 

and the system of equation 

x(t ) = Ax + Bu becomes (8) 

x{t) = Ax + B(-kx ) . (9) 



Lyapunov Tuned State Controlled PWL Model w/ Pl=-2.1 and P2=-l 
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Figure 9 - Lyapunov Optimized Controller 



RICCATI OPTIMIZED CONTROLLER 



Another well-known approach to determining the 
state-variable feedback (SVFB) control gain k is the 
Riccati method to optimize the linear quadratic 
equation (LQE) by minimizing the cost function 
[24,25]: 

/ = Jo ° 0 [* r 7 x + ±u 2 ]dt. ( 12 ) 

This method differs from the Lyapunov method by 
not manually selecting poles by selecting tuning 
weight parameters Q, for system stability, and R, 
sometimes referred to as ^ , for input control effort. 

Using these parameters and minimizing the cost 
function minimize the error thus creating a stable 
system. For simplicity, Q can be set to the I matrix of 
the appropriate size and R, or y, adjusted to give the 
desired results [21,25]. Substituting the input from 
equation (7) yields: 

) = ^"[x T {Q + K T - y K)x\dt. (13) 



This equation is used to optimize the controller, thus 
minimizing the error and ensure stability. Suppose 
there exists a matrix P that satisfies the following 
equity: 



£(x t Px) = -x T (Q + K T ± K)x. (14) 



Substituting of (14) into (13) yields: 
J = * T (0)P*(0) 



(15) 
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where the closed loop system’s stability is assured by 
P remaining constant and x(t) going to zero as t goes 
to infinity. By differentiating (14) substitution into 
(9) and algebraic manipulation, the algebraic Riccati 
equation (ARE) is found: 

A t P + PA + Q- 2 yPBB T P = 0 . (14) 

By selecting y and Q is positive definite matrices and 
solving for the P matrix, the system is guaranteed 
asymptotic stable in its closed loop operation with a 
gain of K = 2 yB T P. These results are shown in 
figure 10. 



Riccati Tuned State Controlled PWL Model w / Q=l and R= 001 
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4. RESULTS 

Figure 11 shows the results of the nonlinear model 
with PID control (a), the PWL model with PID (b), 
the PWL model with Lyapunov tuned poles at -1 and 
-2 (c) and the PWL model with Riccati tuned poles 
with R=0.001 (d). As shown in Figure 3, (a) and (b) 
are very similar, with only a few RPMs of error 
across the operating range and a settling time of about 
5 seconds. Run (c) with the Lyapunov tuned poles 
controller reduces the settling time to about 2 
seconds, but there is a significant amount of 
overshoot during the first second, which causes the 
clutch rapidly lockup and then release, causing a 
significant amount of driveline shock. The Riccati 
tuned controller in run (d) showed the best results 
with near instantaneous response to the control input 
and less than % % of overshoot. This superior 
performance comes at a cost. The closed loop 
feedback gains are in the 10,000s. Depending upon 
the application and hardware, this may not be 
possible, resulting in a slower response with the gains 
in the 10s and 100s. 



ROBUSTNESS 

Another area evaluated in this paper is the 
controller’s robustness to noise and changing inputs. 
The various controllers’ ability to track are shown in 
figure 12 (25%) and figure 13 (45%) with noise at % 
Hz added to the measurement of the output speed 
sensors. It is clearly seen in all figures that the PID 
controller is too slow to react to the RMS of the 
added noise but rather compensates for the minimum 
values of the combined inputs. This anomaly 
exaggerates the error, causing poor system 

performance, while the Lyapunov and Riccati 
controllers are able to react to the RMS of the signal 
thus minimizing the effect on the system’s output. 
Figure 10 (c) shows the system has reached the limit 
where switch gain controllers can operate and they 
become unstable at these high noise levels 



Figure 10 - Riccati Optimized Controller 
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Figure 1 1 - TCC Input and Output Speeds Response to 100 RPM 
Delta Target: 

a) PID Controlled Nonlinear Model, b) PID Controlled PWL 
Model, c) Lyapunov Optimized PWL Model, 
and d) Riccati Optimized PWL Model 
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Figure 12 - TCC Input and Output Speeds Response to 100 RPM 
Delta Target with 25% Measurement Noise: 
a) PID Controlled Nonlinear Model, b) PID Controlled PWL Model, 
c) Lyapunov Optimized PWL Model, and d) Riccati Optimized PWL 
Model 
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1446 


3 


0.751 


4.5 


1505 


1446 


4 


0.752 


4.5 


1505 


1446 


4 


0.753 


4.5 


1505 


1446 


4 


0.754 


4.5 


1504 


1446 


4 


0.755 


4.5 


1503 


1446 


4 


0.76 


4.5 


1502 


1446 


5 


0.765 


4.5 


1499 


1446 


5 


0.77 


4.5 


1496 


1446 


5 


0.775 


4.5 


1491 


1446 


5 


0.78 


4.5 


1467 


1445 


5 


0.7825 


4.5 


1457 


1444 


6 


0.785 


4.5 


1450 


1444 


6 


0.7875 


4.5 


1445 


1443 


6 


0.788 


4.5 


1445 


1443 


6 


0.7885 


4.5 


1444 


1443 


6 


0.789 


4.5 


1443 


1443 


7 


0.8 


4.5 


1443 


1443 


7 


0.85 


4.5 


1443 


1443 


7 


0.9 


4.5 


1443 


1443 


7 


1 


4.5 


1443 


1443 


7 



Table 1 - Nonlinear Model Response 
and PWL segments 



5. CONCLUSIONS 

Our results show that high order nonlinear system can 
be modeled using Piecewise Linearization with high 
degrees of accuracy and robustness. We also found 
that State controllers require fewer segments to 
accurately model and to control a PWL model in part 
because the PWL model segments where each tuned 
there specific operation regions and also because of 
the way the segments where selected caused general 
lineamess between them. The advantage of the 
Riccati tuned controller over typical Lyapunov pole 
place controllers was shown by fast reactions to this 
high order nonlinear system. In practice, it may not 
be possible to duplicate the full benefit of this 
controller in a TCC because it requires higher gains 
in the 10,000s to overcome the 5 order electro- 
hydraulic control circuit in addition to the nonlinear 
clutch action. The rapid response of the Riccati 
controller may also be undesirable because its near 
infinite slope will create driveline shocks similar to 
the original quantitative state lockup TCC. 

A separate area of mixed results is the controllers’ 
robustness to noise. State optimized controllers do an 
excellent job responding to noise when the signal 
bounds keep the controller within the same operating 
segment, however significant state errors can be 
generated when the controller output or model when 
they are required to rapidly oscillate between 
segments to react to the noisy signal. 
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Abstract 

From the past studies the breakdown strength of SF 6 is 
exceptionally subjected to nonlinear electric fields in 
conjunction with metallic particles, consequently this 
effects provokes drop in the breakdown strength beyond 
its limits. In addition to this the SF 6 gas is a hothouse gas 
rouse the global warming, employing this gas is 
extremely undesirable. As a consequence there is a need 
to find alternate solution to this gas one such alternative 
is to use gas mixtures such as SF 6 /N 2 , SF^Air, SF 6 /C0 2 , 
SF 6 /Kr, SF 6 /Ar having supreme dielectric characteristics. 
In this paper Metallic Particles movement is analyzed, 
compared and presented the results of SF 6 /N 2 , 
SF 6 /C0 2 and SF 6 /Air gas mixtures in a Single Phase Gas 
Insulated Substation. 

Keywords: Gas Insulated Substation, gas mixture, 
metallic particle, hothouse gas. 

1. Introduction 

In the early stages of nineteen seventies the gas insulated 
substation were make available to the end users to get the 
advantage of its size and inconsiderateness to 
atmospheric circumstance. The gas insulated substation 
has outstanding reliability and had given solutions to 
worlds power conveyance challenges [7].SF 6 gas is 
commonly used in the high-voltage electrical engineering 
for the reason that it is a praiseworthy dielectric and has a 
tremendous electrical arc interrupting characteristic, 
making it the most suitable gas for substation 
applications. However, there is just one difficulty; it is a 
hothouse gas with a lofty Global Warming Potential 
GWP [1]. 

Diluting SF6 gas mixtures are of specific interest as 
assured reduction of the environmental effects can be 
achieved. Various characteristic of the insulation 
behavior of SF 6 and N 2 mixtures have been published 
[6]. So much of research has been done on sf 6 gas 
mixtures for replacing the sf6 ,some of those mixture gas 
had given better dielectric characteristics the sf6 but they 
are toxic in nature . Examinations were done on SF 6 /N 2 
gas mixture with lower proportion of SF 6 and higher 
proportion of N 2 . The most favorable process was needed 



to find the optimum solution for fixing the percentage of 
the SF 6 .some author had stated that with 12-25% of SF 6 
in a gas mixture can give same insulating capabilities of 
75-85 % of pure SF 6 . 

It is identified that the dielectric strength of the gas 
insulated substation is affected by the presence of 
metallic particles which are free to move due to applied 
electric field in the compressed gas. These particles may 
have insulating or conducting charterstics which are 
originated during the transportation, assembling of 
substation or during switching actions. Movement of the 
insulating particles does not affect the dielectric strength 
of the gas but the movement of conducting particles leads 
to micro discharges and electric stress of high magnitude 
[3-6] which leads to breakdown of the system. 

The work stated in this paper deals the dynamic of 
behaviour of metallic particles in a gas insulated 
substation with various SF 6 gas mixtures. A software 
program is developed for computing the conducting 
particle movement for different proportions of gas 
mixture at a specific electric field using an advanced C++ 
language. In a single phase gas insulated substation a 
fixed geometry of copper and aluminium particles of wire 
type were considered to determine the radial movement 
of the particle. To solving the motion equation analytical 
method has been adopted for various gas mixtures. 

2. Modelling Technique 

The figure 1 shows the parallel conductor bus of gas 
insulated substation enclosed in an enclosure packed with 
mixture of SF 6 and other gas is considered for the study. 




Figure 1 
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A conducting particle of wire type is supposed to be at 
rest on the surface of the enclosure till a sufficient 
voltage is applied to pick up the conducting particle. 
After gaining sufficient charge and surmounting its own 
weight and drag forces, conducting particle winches up 
and progress in the direction of the applied field. 
Numerous of parameters were considered to determine 
the particle motion equation , for instance restitution 
coefficient , Reynolds number, drag coefficient, viscosity 
of gas and so on .though several authors have expression 
for wire like particles but the equations are taken from 
the work of Felici [2] . 

When a conducting particle is subjected to the electric 
field in a gas insulated substation it attains an electric 
charge and undergoes an electrostatic force. If this 
force surpasses the gravitational force, the conducting 
particle is picked up from the enclosure and a zig zag 
movement is detected[8]. 

Ignoring the image charges due to the existence of 
enclosure grounding, the electrostatic force F e is 
proportional to the charge 

F e=Qw E f 1 



For a wire like particle resting on the surface of the 
enclosure [5, 6, 7], for 1 » r, the charge Q w is 



Q w 




2 



A conducting free particle moving in the applied electric 
field is subjected to Electrostatic force (F e ), force due to 
dra g(F d ) and force due to gravitation of earth (F g ) . The 
free conducting particle motion of mass (m) is computed 
using the equation of motion 




-F„ 



— F„ 



3 



The Reynolds number and the viscosity of the two gas 
mixture is required for calculating the motion equation 
The above equation (3) leads to non-linear equation 
which is solved using iterative technique such as Runge- 
Kutta fourth Order Method. 



3. Simulation of Particle motion 

Analytical method is used for calculating electric field 
and charge acquired using equation (1.1) and (1.2) based 
on work of [2, 5 and 8]. The free conducting particle 
movement given in equation (1.3) is simulated using the 
analytical method. Using advanced C++ programming, 
simulation of motion equation of free conducting wire 
particles with 15mm length and 0.5mm radius is carried 
out in gas insulated substation for 132 kV, 145kV and 
220 kV applied voltages. The enclosure radius of 88 mm 
and conductor radius of 33 mm is considered for the 
simulation. It is assumed that Aluminum and Copper 
particles of wire type were present in the gas insulated 
substation 



SF 6 /C0 2 and SF 6 /Air mixture for 132 KV, 145KV and 
220KV can be known. The maximum movement 
obtained for A1 particles and Cu particles for 132 kV is 
11.876097, 11.627509, 11.889343 and 2.220426, 

2.227485, 2.250473. The maximum movement obtained 
for A1 particles and Cu particles for 145 kV is 14.608354, 
15.656908, 14.845436 and 2.778246, 2.778473, 

2.832792. The maximum movement obtained for A1 
particles and Cu particles for 220 kV is 29.648418, 
30.830299, 30.822065 and 9.429732, 9.444219 and 
9.558602. 



S 40 
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|| 2 . 
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□ SF6/N2 

□ SF6/C02 

□ SF6/Air 



132 KV 145 KV 220 KV 
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Figure 2 




Figure 3 



From the figures 4 and 5 the maximum movement of 
Aluminum and Copper particles for 40% SF6 in SF 6 /N 2 , 
SF 6 /C0 2 and SF 6 /Air mixture for 132 KV, 145KV and 
220KV can be known. The maximum movement 
obtained for Al particles and Cu particles for 132 kV is 
12.087142, 11.644285, 12.737603 and 2.219076, 

2.287596, 2.209517. The maximum movement obtained 
for Al particles and Cu particles for 145 kV is 15.698053, 
14.822464, 14.63584and 2.843645, 2.783404, 2.788052. 
The maximum movement obtained for Al particles and 
Cu particles for 220 kV is 30.600975, 30.605444, 
30.880322 and 9.540246, 9.581879and 5.644269 
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4. Results and Discussion 

From the figures 2 and 3 the maximum movement of 
Aluminum and Copper particles for 20% SF6 in SF 6 /N 2 , 
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Figure 5 



From the figures 6 and 7 the maximum movement of 
Aluminum and Copper particles for 60% SF6 in SF 6 /N 2 , 
SF 6 /C0 2 and SF 6 /Air mixture for 132 KV, 145KV and 
220KV can be known. The maximum movement 
obtained for A1 particles and Cu particles for 132 kV is 
11.897035, 11.897187, 11.892382 and 2.224117, 

2.224119, 2.220146. The maximum movement obtained 
for A1 particles and Cu particles for 145 kV is 15.126985, 
14.842751, 14.844956 and 2.833687, 2.806137, 

2.778574. The maximum movement obtained for A1 
particles and Cu particles for 220 kV is 31.564141, 
31.100422, 31.222407 and 9.631787, 9.473286and 
9.354718 





Figure 7 



From the figures 8 and 9 the maximum movement of 
Aluminum and Copper particles for 80% SF6 in SF 6 /N 2 , 
SF 6 /C0 2 and SF 6 /Air mixture for 132 KV, 145KV and 
220KV can be known. The maximum movement 
obtained for A1 particles and Cu particles for 132 kV is 
11.935179, 12.140633, 11.979409and 2.22292, 2.23103, 
2.288991. The maximum movement obtained for A1 
particles and Cu particles for 145 kV is 14.939284, 
14.873562, 14.988114 and 2.837505, 2.831574, 

2.841872. The maximum movement obtained for A1 
particles and Cu particles for 220 kV is 30.314663, 



30.666139, 31.39572 and 9.552808, 8.310761 and 
9.587192. 




Figure 8 




Figure 9 



From the figures 10 and 11 the maximum movement of 
Aluminum and Copper particles for 100% SF6 in SF 6 /N 2 , 
SF 6 /C0 2 and SF 6 /Air mixture for 132 KV, 145KV and 
220KV can be known. The maximum movement 
obtained for A1 particles and Cu particles for 132 kV is 
11.95468 and 2.192154. The maximum movement 
obtained for A1 particles and Cu particles for 145 kV is 
14.939284 and 2.834764. The maximum movement 
obtained for A1 particles and Cu particles for 220 kV is 
31.384289 and 9.595386. 




Figure 10 




Figure 1 1 
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5. Conclusion 

It is observed that, as the applied voltage increases the 
particle movement increases. It is also observed that the 
Aluminum and Copper particles movement is almost 
identical for 20%, 40%, 60%, 80% SF6 and 100% SF 6 in 
gas mixtures mentioned. In order to reduce the 
environmental effects occurred due to SF 6 gas utilization 
in electrical industry as insulating medium particularly in 
Gas insulated substations, 20% SF 6 in gas mixture can be 
used as insulating medium in the gas insulated substation. 
Among the gas mixtures given above SF 6 /N 2 gives the 
better results in terms of economy, dielectric capability 
and environment issues. 
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el: Particle movement ] 


? or 0% sf6 in gas mixture for different voltages 


S. 

No 
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132 KY 
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movement 
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e 2: Particle Movement for 20% sf6 in gas mixture for different voltages 



s. 

No 


20% SF 6 


132 KV 


145 KV 


220KV 


AL 


CU 


AL 


CU 


AL 


CU 
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Particle 
movement 
in mm 


SF 6 /N2 


11.876097 


2.220426 


14.608354 


2.778246 


29.648418 


9.429732 


2 


sf 6 /co 2 


11.627509 


2.227485 


15.656908 


2.778473 


30.830299 


9.444219 


3 


SF 6 /Air 


11.889343 


2.250473 


14.845436 


2.832792 


30.822065 


9.558602 
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e 3: Particle movement for 40% sf6 in gas mixture for different voltages 



S. 

No 


40% SF 6 


132 KV 


145 KV 


220KV 


AL 


CU 


AL 


CU 


AL 


CU 


1 


Particle 
movement 
in mm 


SF 6 /N2 


12.087142 


2.219076 


15.698053 


2.843645 


30.600975 


9.540246 


2 


sf 6 /co 2 


11.644285 


2.287596 


14.822464 


2.783404 


30.605444 


9.581879 


3 


SF 6 /Air 


12.737603 


2.209517 


14.63584 


2.788052 


30.880322 


5.644269 
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A1 Particle,40% SF6 : 



Cu Particle,40% SF6 : 





132 KV 145 KV 



220KV 



□ SF6/N2 
■ SF6/C02 

□ SF6/Air 



Tab: 



e 4: Particle movement for 60% sf6 in gas mixture for different voltages 



s. 

No 


60% SF 6 


132 KV 


145 KV 


220KV 


AL 


CU 


AL 


CU 


AL 


CU 


1 


Particle 
movement 
in mm 


SF 6 /N2 


11.897035 


2.224117 


15.126985 


2.833687 


31.564141 


9.631787 


2 


sf 6 /co 2 


11.897187 


2.224119 


14.842751 


2.806137 


31.100422 


9.473286 


3 


SF 6 /Air 


11.892382 


2.220146 


14.844956 


2.778574 


31.222407 


9.354718 



A1 Particle,60% SF6 : 




Cu Particle,60% SF6 : 



□ SF6/N2 
■ SF6/C02 

□ SF6/Air 




□ SF6/N2 
■ SF6/C02 

□ SF6/Air 



Tab 



e 5: Particle movement for 80% sf6 in gas mixture for different voltages 



S. 

No 


80% SF 6 


132 KV 


145 KV 


220KV 


AL 


CU 


AL 


CU 


AL 


CU 


1 


Particle 
movement 
in mm 


SF 6 /N2 


11.935179 


2.22292 


14.939284 


2.837505 


30.314663 


9.552808 


2 


sf 6 /co 2 


12.140633 


2.23103 


14.873562 


2.831574 


30.666139 


8.310761 


3 


SF 6 /Air 


11.979409 


2.288991 


14.988114 


2.841872 


31.39572 


9.587192 



A1 Particle.80% SF6 : 




Cu Particle.80% SF6 : 



□ SF6/N2 
■ SF6/C02 

□ SF6/Air 




□ SF6/N2 
■ SF6/C02 

□ SF6/Air 
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Table 6: Particle movement for 100% sf6 in gas mixture for different voltages 



s. 

No 


100% sf 6 


132 KY 


145 KY 


220KV 


AL 


CU 


AL 


CU 


AL 


CU 


1 


Particle 
movement 
in mm 


SF 6 /N2 


11.95468 


2.192154 


14.695567 


2.834764 


31.384289 


9.595386 


2 


sf 6 /co 2 


11.95468 


2.192154 


14.695567 


2.834764 


31.384289 


9.595386 


3 


SF 6 /Air 


11.95468 


2.192154 


14.695567 


2.834764 


31.384289 


9.595386 



A1 Particle. 100% SF6 : 




Cu Particle. 100% SF6 : 



□ SF6/N2 
■ SF6/C02 

□ SF6/Air 




□ SF6/N2 
■ SF6/C02 

□ SF6/Air 
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Abstract 

In this work, a simplified technique to control and 
stabilize the passively Q-switched laser pulses is 
proposed and denoted by automatic pulsed pumping 
technique. The technique is applied to a diode-pumped 
Yb:YAG laser passively Q-switched by Cr 4+ :YAG as 
saturable absorber by numerically solving the coupled 
rate equations describing the system. The technique is 
based on the detection of the generated Q-switched 
pulses which alter the pumping power either to be 
switched off or set to a lower value depending on the 
selected modulation depth. 

The system was pumped using different pumping 
techniques; continuous, pulsed pumping and automatic 
pulsed pumping in order to compare the time stability of 
the generated Q-switched pulses. The proposed technique 
has shown a good control of the pulse repetition rate of 
Q-switched pulses with high stability. The technique 
allowed high pulse repetition rate as similar to the pulsed 
pumping scheme while maintained the time stability of 
the continuous pumping scheme. 1 

Keywords: Q-switched lasers; pulse repetition rate; time 
stability; pulsed pumping. 

Nomenclature: 

cp the photon density in the laser cavity, 

t r the cavity round-trip time, 

o the stimulated emission cross-section of the 

laser crystal, 

o g the absorption cross-section of ground state of 
the saturable absorber, 

G e the absorption cross-section of the excited state, 
N the population inversion density of the gain 

medium, 

n g and the absorber ground state and excite state 

n e population density respectively, 

n s0 the total population density of the saturable 



1 This study has been implemented at National Institute 
of Laser Enhanced Sciences, Cairo University, Giza, 
Egypt. 



absorber, 

/ the length of the laser crystal, 

l s the length of the saturable absorber, 

l! the optical length of the laser cavity, 

R the reflectivity of the output coupler, 

L nonsaturable intracavity round-trip dissipative 

optical loss, 

W p the volumetric pump power, 

c the speed of the light, 

y the inversion reduction factor, 

r the lifetime of the upper laser level of the gain 

medium 

r s excited-state lifetime of the saturable absorber, 

Wp n (t) the normalized volumetric pump power , 

(j) n (t) the normalized photon density, 

N n (t) the normalized population inversion density of 
gain medium, 

md the modulation depth of the pulsed pumping 
power, 

T the period of pulsed pumping power. 

1. Introduction 

Diode pumped Q-switched lasers are currently used to 
produce pulses with high peak power at high repetition 
rates and short nanosecond pulse widths. These types of 
lasers are used in different applications such as medicine, 
material processing and military. 

Q-switching can be established by either active or passive 
techniques. Although the active Q-switching provides 
higher power extraction, higher control on the repetition 
rate and lower time jitter, it is more complex and 
expensive than the passive Q-switching [1]. 

In the past, passive Q-switch element was first 
implemented by organic dye cells [2]. Due to poor 
durability of dye cell (degradation of dye) and low 
thermal limit of plastic material, it was replaced by doped 
crystals which provide higher power and better reliability 
[3]. The first new material was F 2 :LiF color center crystal 
[2]. Then, Cr-doped crystals showed more photochemical 
and thermal stability [4]. Different crystals were involved 
with Chromium ions (Cr 4+ ) such as Cr 4+ :YAG, 
Cr 4+ :GSGG and Cr 4+ :YSO as they have favored 
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characteristics such as large absorption cross section, low 
saturable intensity at the lasing wavelength. Therefore, 
the Cr-doped crystals became promising Q-switching 
elements due to their low cost, long life, compactness. 
Cr 4+ :YAG is the most of attractive saturable absorber due 
to the favorable chemical, thermal and mechanical 
properties of YAG host as well as its high absorption 
cross section of laser wavelengths centered at 410, 480, 
640 and 1050 nm [5-6]. 

Other solid-state saturable absorbers such as 
Semiconductor Saturable Absorber (SESAM), GaAs 
wafers, and Saturable Bragg Reflector (SBR), are also 
available and investigated in many previous works [7-8]. 
Thence, solid-state Q-switched lasers are a good 
candidate in many applications due to their compactness, 
flexibility and reliability. 

Most of previous work have studied diode pumped Nd- 
based lasers using Cr 4+ :YAG as a Q-switch element. This 
passive Q-switching combination was very attractive due 
to their highly stable phototropic centers and thermal 
conductivity [9]. Lot of work have been reported for 
modeling of passively Q-switched lasers with Cr 4+ :YAG 
as saturable absorber to determine pulse energy, peak 
power, pulse width, polarization, and thermal behavior 
[4,10-14]. 

The stabilization of pulse to pulse time, also referred to 
as time jitter, is a crucial requirement for many 
applications such as remote sensors, range finders or 
three dimensional scanners, Laser Induced Breakdown 
Spectroscopy (LIBS), particle flow measurements and 
defense [3,5,9,15-18]. Therefore, different techniques 
have been proposed to minimize the time jitter. These 
techniques were either introducing different operating 
conditions for active medium, saturable absorber and 
cavity length or modifying the system configuration. 
These techniques will be discussed in more details in 
Section 2. 

In the last decade, passively Q-switched Yb:YAG lasers 
using Cr 4+ :YAG have been reported due to its simple 
energy level scheme, long upper laser level lifetime 
(951psec) [7,19]hich is suitable for storing energy and Q- 
switched pulse generation. In addition, it has low 
quantum defects which decrease heat generation during 
lasing three times less than Nd-based laser systems [20] 
and hence higher pumping efficiency [5]. The lower 
stimulated emission cross section of Yb:YAG increases 
the saturation fluence and the stored energy required for 
bleaching the Q-switch element which yield to shorter 
pulse length higher energy pulse and higher extraction 
efficiency. Besides, the large absorption cross-section 
tolerates the pumping laser wavelength [5,21]. Q- 
switched Yb:YAG laser was also realized by using GaAs 
wafer as saturable absorber [7,22]. 

In this work, the Yb:YAG laser Q-switched by Cr 4+ :YAG 
as saturable absorber was numerically investigated. The 
coupled rate equations describing the photon density, 
population inversion density of the active medium and 
the ground state population density of the saturable 
absorber were solved using 4 th order Runge-Kutta 
algorithm. Using both continuous and pulsed pumping 
(PP) techniques, the PRR and stability were determined. 
To control the PRR and maintain high stability of the Q- 



switched pulses, an enhanced pulsed pumping technique 
was proposed to automatically control the pump power. 

In section 2, the literature review of the studies 
discussing the different techniques to reduce the time 
jitter are presented. The rate equations of the Yb:YAG 
laser Q-switched by Cr 4+ :YAG as saturable absorber are 
provided in section 3. The results of continuous and 
pulsed pumping are discussed in Section 4. Then, the 
proposed technique of automatic pulsed pumping is 
represented as well as the resulting PRR and time jitter. 

2. Literature Review 

As it was mentioned in Section 1 , many studies discussed 
the control of pulse repetition and time jitter through 
different techniques and schemes. 

Many factors contribute to the pulse to pulse time jitter 
such as the environmental and pump instabilities 
including the intensity variations, changes in pump 
spatial and spectral properties as well as the spontaneous 
nature of passive Q-switching [1] [23]. The 

environmental and pump instabilities are contingent on 
the stability of the surroundings and the reliable pumping 
electronics and laser diodes. Environmental isolation can 
reduce the time jitter by factor of 5 in the range of psec 
[24], 

The spontaneous noise is also a major source of 
instability since the pulse formation depends on the 
spontaneous emitted photons when the photon flux is 
close to the lasing condition [3] [23]. Lor a stable 
pumping, the stimulated emission is more promoted over 
the spontaneous one. Therefore, for higher pumping 
power, the pulse generation can be more stabilized. 

Lai et al. had studied the effect of slightly modulating the 
pump diode power [15]. The repetition rate of a Nd:YAG 
laser passively Q- switched by Cr 4+ :YAG was 
experimentally shown to reach a time stability of 10" 6 
over 10 6 pulses, short-time period jitter about 1% in less 
than 1 sec, and a slow variation of 6 % over a few minutes 
due to environmental thermal drifts. 

However, the modulation of pump power may result in 
some undesirable dynamical behavior such as periodic 
and chaotic oscillations [25]. The study carried out by K- 
G Hong et al, was based on Nd:YV0 4 laser passively Q- 
switched by Cr 4+ :YAG. The stabilization of pulse 
repetition rate and time jitter was controlled by careful 
selection of the modulation frequency to match the 
repetition rate of free running mode otherwise periodic 
and chaotic behavior may occur. 

Also, the effect of modulating the pumping laser diode 
was also studied in microchip structured Yb:YAG Q- 
switched laser as well as introducing a Brewster polarizer 
in order to maintain linear polarization [5]. The time 
jitter was drastically decreased from psec to nsec range. 
An alternative technique to control the pulse repetition 
rate and time jitter was changing the transmission of the 
saturable absorber by optical bleaching. Klimov et al. [9] 
experiment were conducted using a diode pumped 
Nd:YAG laser while controlling the orientation of the Q- 
switch element around its optical axis. 

Another technique of saturable absorber bleaching is 
applying external optical modulation to the saturable 
absorber [3]. The modulated output of 910 nm laser diode 
was injected into the Cr:YAG saturable absorber to alter 
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its population difference and hence its transmission. The 
modulation frequency has been chosen to match the free 
running repetition rate. The time jitter was reduced by a 
factor of 2.5. 

Instead of applying optical modulation, a direct bleaching 
of the Cr:YAG saturable absorber was involved [1], A 
1036 nm laser diode bar was driven to trigger a 920 W 
bleaching pulse of 15 psec to turn the cavity loss of the 
saturable absorber to a lower lasing threshold. With 
careful selection of the triggering time, the time jitter was 
reduced by factor of 12 from 241 nsec to 20 nsec. 
Moreover, an active Q- switched laser bleaching the 
Cr:YAG saturable absorber of a passively Nd:YAG Q- 
switched laser has been reported [26]. Sub-nanopulses 
with 632 psec time jitter have been realized. 

The previously mentioned techniques were based on 
either modulating the pumping laser diode or optical 
bleaching of the saturable absorber. S. Su [17] had 
proposed to merge the two techniques while considering 
the fast switching of the laser diode. Operating the laser 
diode in two conditions; less than and close to the 
threshold value as a prepumping of the active medium, 
the time jitter is reduced from the range of psec to 150 
nsec. 

The reflectivity and the length of the lasing cavity are of 
the design parameters that can contribute to the control of 
the Q-switched output pulse [23]. Kausas et al, have 
proposed to extend the lasing cavity with a coupled 
cavity of the same length of the initial cavity to match 
their free spectral range. The coupled cavity is realized 
by an additional mirror mounted on a piezoelectric 
transducer. By changing the position of the additional 
mirror, to switch the gain thresholds between high and 
low, the pulse generation can be controlled and hence the 
time jitter (from 190 nsec to 24 nsec). 

The active medium of Q-switched lasers can play a role 
in the control of the output pulses. Some authors changed 
the host of active ions from Nd:YAG to Nd:YV0 4 [17] 
and from crystalline form to ceramic one [26,27]. In this 
work, Yb:YAG active lasing medium is chosen. 

For further control of the pulse generation, the systems 
can migrate to more complicated configurations such as 
using master/slave lasing cavities [28], hybrid active and 
passive Q-switched lasers [17], mixing Q-switched with 
mode-locked mechanisms [29] and moving to optical 
fiber laser configurations [30]. 



3. Rate Equations 

According to the passively Q-switched laser theory 
[4,10,31], the modified coupled rate equations of photon 
density and population inversion density of gain medium 
and the saturable absorber in the passively Q-switched 
resonator, which include the excited-state absorption of 
the saturable absorber, the pump term and the population 
reduction factor of the laser, are given by: 



dcp 

dt 



2aNl — 2 OgUgls ■ 



-In 



Q 



dN n 

_ = _ yaC(pN __ +Wp 



( 1 ) 

( 2 ) 



dn, 

dt 



■g n sO n g 

= —OgCCptlg 



( 3 ) 



n g + n e = n s0 (4) 

where cp is the photon density in the laser cavity, 
t r the cavity round-trip time (t r =2 tic ) of laser 

cavity of optical length t , 

o the stimulated emission cross-section of the 

laser crystal, 

o g the absorption cross-section of ground state 

of the saturable absorber, 

o e the absorption cross-section of the excited 

state, 

N the population inversion density of the gain 

medium, 

n g andn e the absorber ground state and excite state 
population density respectively, 
n s0 the total population density of the saturable 

absorber, 

/ the length of the laser crystal, 

l s the length of the saturable absorber, 

R the reflectivity of the output coupler, 

L the nonsaturable intracavity round-trip 

dissipative optical loss, 

W p the volumetric pump power, 

c the speed of the light, 

y the inversion reduction factor ( y =2 for Yb 3+ - 

doped three-level solid state lasers), 
r the lifetime of the upper laser level of the 

gain medium 

r s the excited-state lifetime of the saturable 

absorber. 



The relationship between the population density in the 
atomic level ‘ 1’ (n s i) and the total population density (n s0 ) 
of the saturable absorber can be given by [4] 

n si/ n so = 0V / UiT (5) 

where a describes the relative strength of the absorption 
transitions of the laser active medium and the saturable 
absorber ( a=<j g ly<j ) and n t is the threshold population 
inversion density in the active medium at which the Q- 
switched pulse is initiated, given by 

\n(l/R)+2o g n s0 l s +L 

n - = d (6) 

By solving the rate equations altogether, the population 
inversion density ( n t ) at peak power and the final 
population inversion density (nj) of the obtained Q- 
switched pulse can be described by the following 
transcendental equations: 

n £ = n £ o + / 1 _n £ oU^“ 

rii rii \ rti ) \riiJ 

rij rij Vrij/V rij/aV Vn^/ / (8) 

where 8= cr e /cr g and n t0 is defined by 



n t o = 



\n(l/R)+S ln(l/To )+L 



2al 



The peak power P(t ) then can also be determined by 

P ^) = ^7 ln (fi) Vmax 
While the laser pulse energy ( E ) is expressed as 



( 9 ) 



(10) 
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Thus, the pulse width can be calculated from 

t p = E/P(t ) (12) 

By solving Equation (1) at steady-state, the total loss can 
be given by: 

2 a g n g l s + 2 a e n e l s + ln(l /R) + L 

Loss = - — - I 1 -*) 

2al 

In the case of continuous pumping, the pumping term W P 
is considered as a constant value. While in the case of 
pulsed pumping, the pumping term will be modified to a 
rectangular pulse train as shown in Figure 1, where T is 
the period of pumping cycle, p w is on duty time and md is 
the modulation depth. However, the continuous pumping 
can be regarded as a special case of pulsed pumping 
when md= 0. 

The time between two successive Q-switched pulses can 
be denoted by t pp . For N pulses, the PRR can be 
calculated from the average of these times x av as 
1 

PRR = — (14) 

^av 

The time jitter between pulses can be also obtained by 
calculating the time deviation of the generated Q- 
switched pulsed from the average time x av . The time jitter 
can be considered as an indication of the time stability of 
the generation of the Q-switched pulses. 

jitter = <|x pp — x av |) (15) 

Hence, the time instability can be determined by the ratio 
between the jitter and the average time, 

instability = ^ Tpp Tav ^ (16) 

T av 




Figure 1. Rectangular Pulse train form of pulsed pumping in PP 
technique 



4. Simulation Results 

The coupled rate equations (1) to (3) are solved 
numerically using 4 th order Runge-Kutta method. The 
spectral parameters of 10% at Yb:YAG passively Q- 
switched using Cr 4+ :YAG are listed in Table (1) [32]. 
Figure 2 shows the numerical stimulation results of the 
first four pulses. The first pulse was obtained after 197 
psec, while the following pulses were generated with 
average time T av of 202 psec and the free running pulse 
repetition rate of 4.95 KHz. The peak power can be 
determined from Equation (12), where the pulse width is 
15.2 ns (from Figure 3) and pulse energy is 117 pJ (from 
Equation (11)). Thus, the calculated peak power is 7.7 
KW. 

Figure 3 illustrates the first generated Q-switched laser 
pulse. When the photon density inside the laser cavity 



increases, the loss decreases accordingly as a result of the 
bleaching effect of Cr 4+ :YAG saturable absorber. When 
the population inversion density becomes equal to the 
cavity loss, the photon density attains its peak value then 
decreases as long as the population inversion density is 
decreased. Figure 4 represents the normalized volumetric 
pump power (W Pn (t)), the normalized photon density 
(<j) n (t)) and the normalized population inversion density of 
gain medium (N n (t)) in the case of md= 0, which 
corresponds to the continuous pumping case. The 
calculated PRR is ~ 5 kHz while the jitter is 3.2xl0" n sec 
corresponding to time instability of 1.6x1 0" 7 . 



Table 1. The spectral parameters of Yb:YAG passively Q-switched 
using Cr 4+ :YAG . 



Symbol 


Value 


<7 


2.3xlO -20 cm 2 @ 1030 nm 


/ 


0.4 cm 


c 


3xl0 10 cm/sec 




3.2xl0" 18 cm 2 @ 1.06 pm 


Oe 


4.5xl0" 19 cm 2 @ 1.06 pm 


Is 


0.94 mm 


n s o 


0.87xl0 17 ions/cm 3 


R 


85% 


To 


97% 


V 


5.2 cm 


y 


2 


K 


1.12xl0 23 sec" 1 cm" 3 


ht)p 


2.1 12x1 0" 19 Joule 


T 


951 psec 


T s 


3.4 psec 


hv 


1.93xl0" 19 Joule 


Cyb 


13.8xl0 20 ions/cm 3 (10 at. ( 


L 


8% 


a 


10 cm' 1 



The effect of pulsed pumping on the generation of output 
pulses is studied for different pumping parameters; the 
period T, the duty cycle p w and the modulation depth md. 
Since the average time between pulses is larger than 200 
psec, a pulsed pumping of period T=300 psec and duty 
cycle of 250 psec was first applied. Then, the period was 
fixed and the duty cycle was decreased to 200 psec then 
175 psec. All the different cases are studied for different 
modulation depths md= 0.25, 0.5, 0.75, 1. A special 
condition of modulation depth is when the pumping is 
operated in ON-OFF pumping scheme. 

Figure 5. a shows the ON-OFF pumping case. Since the 
ON-time of the first period is 250 psec which is larger 
than the first Q-switched laser pulse produced in the 
continuous pumping case, the first generated pulse is also 
generated at 202 psec. 

Due to the discontinuity of the pumping (OFF-time) 
during 50 psec, the next pulses are generated with T av 260 
psec, PRR of 3.8 kHz and jitter 1.56xl0" 5 sec 
corresponding to time instability of 6%. Using different 
modulation depths, (m<7=0.25,0.5,0.75 and 1), it can be 
noticed that the PRR has increased up to its double value 
for md= 1 while the stability has shown a degradation as 
the jitter increased by order of 6 with respect to 
continuous pumping 
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Figure 2. Population inversion density N(t), total loss of laser cavity 
(Loss) and photon density (cj>(t)) 




time [jis] 

Figure 3. Population inversion density N(t), total loss of laser cavity 
(Loss) and photon density, (cj>(t)) at the first Q-switched pulse. 
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Figure 4. The normalized volumetric pump power (W Pn (t)), the 
normalized photon density (c|> n (t)) and the normalized population 
inversion density of gain medium (N n (t)) in the case of continuous 
pumping (md= 0) 

The variation of PRR, jitter time and instability can 
be explained from the temporal behavior of population 
inversion density of gain medium N n (t) for each 
modulation depth, as shown in Figure 5. a After 
generation of the first Q-switched pulse, the remaining 
ON-time in the duty cycle allows the growth of the 
population inversion density but cannot acquire the 
threshold value (ni) before the end of ON-time duration. 
During the OFF-time, the population remains almost 
constant until the beginning of the second period of the 
pumping and the population restarts to grow to generate 
the second Q-switched pulse. 

The rate of increase of the population inversion density 
during the OFF-time is inversely proportional to the 
modulation depth; as the modulation depth increases, the 
value of the pumping is decreased. Therefore, the rate at 
which population inversion density increases is slower 
than that of the smaller modulation depths. The case of 




ON-OFF pumping represents a special case of pulsed 
pumping, into which during the OFF-time, the pumping 
is turned off and the population inversion density 
decreases rather than increases as in other pulsed 
pumping cases in which the pumping is not turned off. 
Depending on the modulation depth of the pulsed 
pumping, the time required for the generation of the Q- 
switched pulse decreases with the increase of the 
modulation depth and the pulse repetition rate is 
increased. Therefore, the ON-time of the pumping may 
allow generating more than one Q-switched pulse in the 
same period. 

Keeping the period of pulsed pumping r=300psec, while 
varying the duty cycle from 250 psec to 200 psec then to 
175 psec, the PRR and time jitter have been studied and 
the results are shown in Figures. 6 and 7, respectively. It 
can be seen that PRR is proportional to the modulation 
depth as well as the duty cycle duration. On the other 
hand, the stability decreases drastically for ON-OFF 
pumping cases then becomes more stable in md= 0.25 
case and decreases again when m<7=0. 5,0.75 and 1. 
However, the jitter increases with the decrease of the 
duty cycle time and the stability increases. 

The effect of the period of the pulsed pumping has been 
also studied. The period is decreased from 300 psec to 
150 psec with different duty cycles. The resulting PRR 
and time jitter values for all cases are summarized in 
Table 2. 

It can be seen that within the same pump period, as the 
duty time decreases the PRR decreases due to decreasing 
of the pump time and results in the decrease of the 
amount of optical pumping delivered to the active 
medium during the pumping period, (i.e. decrease the 
average pump power). This occurs for all modulation 
depths. For example if the period is 250psec, and the 
modulation depth is 0.5, the PRR decreases from 5.8 to 
4.94 kHz as the duty time decreased from 200 to 125 
psec. 

Also, it can be observed that within the same pumping 
period and the same duty cycle, the PRR increases with 
increasing the modulation depth since more pumping 
power is gained by the active medium. If 250 psec is 
chosen and 150 psec duty time, the PRR increases 5.1 to 
5.78 kHz when increasing the modulation depth. 

If the duty time is remained constant and the period is 
decreased, the PRR increases since the OFF-time 
decreases and the average pump power delivered to the 
active medium increases. For a 200 psec duty time and 
modulation depth 0.75, the PRR is increased from 5.93 to 
7 kHz as the period is decreased from 300 to 225 psec. 

On the other hand, if the time jitter is examined, it can be 
seen that with the same pump period the time jitter 
increase with decreasing the duty time which means that 
the system is less stable. This behavior is maintained for 
all modulation depths. While if the pump period and duty 
time are kept constant and time jitter increases with the 
increase of modulation depth. Keeping a constant duty 
time, the time jitter is decreasing with the pump period 
decrease for all modulation depths. In general, it can be 
seen that the minimum values of time jitter occur at 
modulation depth 0.25 while the PRR is close to the free 
running case. This can be explained that the change of 
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Figure 7. The time jitter versus modulation depth (T=300psec) 



The pumping power with small values keep the situation 
as close as the original case. 

To summarize, the minimum time jitter was obtained 
with period 200 psec, duty cycle 100 psec and 
modulation depth md= 0.25 corresponding to time 
instability of 5.5 xlO" 4 , but the PRR is less than free 
running value. The highest PRR can be obtained when 
the pulsed pumping has period T = 225 psec, duty cycle 
200 psec and modulation depth md= 1, however, the time 
jitter is increased to 9.27 psec and corresponding 
instability 0.07. 

Although the pulsed pumping technique can be used to 
increase the PRR, however, the degradation of the 
stability was a drawback. In order to increase the stability 
of the output pulses, we proposed a new scheme to 
control the pumping so that once the laser pulse is 
generated, it is detected and the pumping is stopped for a 
previously determined off time. 

After the off time duration, the pumping is restarted to 
attain the required the population inversion density and 
the process of pulse generation is to recommence. This 
technique can be denoted by Automatic Pulsed Pumping 
(APP) technique as shown in Figure 8. Herein, a lOpsec 
OFF-time is proposed and modulation depths vary from 
0.25 to 1 by steps of 0.25. In addition, the ON-OFF 
pumping case is also considered. 

In the first case of ON-OFF pumping, the first Q- 
switched laser pulse is generated at 217psec as shown in 
Figure 9. a. Once the pulse is detected, the pumping is 
switched off for the predetermined time (lOpsec) into 
which the detection and control times are considered. 
After the OFF-time, the pumping is switched ON to 
generate the next Q-switched pulse and the process will 
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Figure 8. The proposed setup of automatic pulsed pumping system 



Table 2. The PRR and time jitter of pulsed pumping technique. 



T (ns) 


Duty 

(ns) 


ON-OFF 


md= 


0.25 


md 


=0.5 


md= 


0.75 


md=l 


PRR 

(kHz) 


jitter 

(10- 6 ) 


PRR 

(kHz) 


jitter 

(10 6 ) 


PRR 

(kHz) 


jitter (10' 
6 ) 


PRR 

(kHz) 


jitter 

(10*) 


PRR 

(kHz) 


jitter 

(10*) 


300 


250 


3.96 


15.6 


5.43 


5.35 


5.9 


10.2 


6.67 


13.1 


7.23 


18.6 




200 


3 


27 


5.2 


8.75 


5.4 


17 


5.93 


26 


6.2 


35 




175 


2.5 


67.5 


5.06 


9.3 


5.2 


18.8 


5.55 


30.3 


5.68 


37.6 


250 


200 


3.78 


9.1 


5.4 


4.2 


5.8 


8.75 


6.58 


13.8 


7 


16.3 




150 


2.65 


59.7 


5.1 


5.9 


5.23 


12.3 


5.6 


19.8 


5.78 


28.6 




125 


2 


2.1 


4.95 


5.9 


4.94 


11.7 


5.13 


18 


5.15 


26.7 


225 


200 


4.3 


3.75 


5.5 


1.91 


6.1 


4.1 


7 


6.7 


7.6 


9.27 




175 


3.58 


21.4 


5.34 


2.97 


5.74 


6.9 


6.46 


12.5 


6.87 


15.5 


200 


150 


3.48 


29.8 


5.3 


1.31 


5.66 


4.1 


6.32 


8.47 


6.67 


12 




100 


2.02 


67 


4.99 


0.11 


5.03 


1.27 


5.01 


2.45 


5.02 


4.9 


150 


100 


2.98 


24.7 


5.2 


3.65 


5.42 


6.2 


5.94 


7.48 


6.21 


5.68 




75 


2.02 


41.1 


4.94 


4.84 


4.94 


9.85 


5.15 


14 


5.15 


18.8 



continue automatically. The PRR and the jitter of the 
generated Q-switched pulses are calculated. The PRR is 
4.4 kHz and the jitter is about 2.1 xl0" n sec. Although the 
PRR has been decreased, it can be seen that the time jitter 
is in the same order of the continuous pumping and the 
time instability is in the order of 10" 7 . 

When the modulation depth of the automatic pulsed 
pumping varies from 0.25 to 1, the PRR is increased to 
7.14 kHz as shown in Figure 10. On the other hand, the 
jitter is slightly affected but is maintained within the 
same order of 10" n sec (Figure 11). Therefore, the 
automatic pulsed pumping has combined the advantages 
of controlling the PRR while keeping the system stability. 
This can be explained as during the ON and OFF timing 
of pulsed pumping, the photon flux and population 
inversion density are not starting and ending at the same 
values. While in the automatic pulsed pumping, the OFF- 
time starts just after the detection of the Q-switched pulse 
and dominates for a predetermined time. Thus, the 
population inversion density exhibits the same decay and 
growth before starting a new pumping period in the case 
of ON-OFF pumping or different modulation depths, 
respectively. Hence, each pulse is subjected to similar 



initial conditions and exhibits the same response 
providing stable Q-switched pulses. 
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5. Conclusion 

In conclusion, the pulse repetition rate PRR, time jitter 
and the time instability of the Q-switched laser pulses 
have been numerically investigated when different 
pumping techniques are applied. The system under 
investigation was Yb:YAG laser passively Q-switched by 
Cr 4+ :YAG as saturable absorber. The free running PRR is 
close to 5 kHz, the peak power 7.7 kW, pulse energy 117 
pj and pulse width 15.2 nsec. 

The PRR can be controlled by using pulsed pumping 
technique. The PRR can be increased from 5 kHz in the 
case of continuous pumping to 7.6 kHz when different 
modulation depths are applied. However, the stability of 
the laser pulses has been decreased since the time jitter 
increased by order of 6 (from 10" 11 to 10" 5 sec). 

A simplified technique has been proposed to control the 
PRR as well as keeping high stability. The technique is 
based on the automatic control of the pulsed pumping so 
that once the Q-switched pulse is generated, the system 
automatically switches the pump status. The pumping 
power is switched off in the case ON-OFF pumping or 
decreased according to the modulation depth which 
varies from 0.25 to 1. Hence, the technique can be 
nominated by automatic pulsed pumping technique. A 
PRR of 7.14 kHz has also been obtained while 
maintaining the time jitter in the order of 10" 11 sec as in 
the continuous pumping case. Therefore, the system 
operates in high repetition rate while keeping the stability 
of pulse to pulse generation. 
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As future work, the numerical model can be extended to 
include the spatial distribution and the longitudinal 
variation of the photon density in order to investigate the 
thermal behavior for proposed pumping techniques and 
compare it with the previous work for continuous 
pumping. In addition, the automatic pulsed pumping 
technique can be experimentally implemented in order to 
validate the results of the numerical simulation. 
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